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INTRODUCTION 


Volume II, Computer Program Documentation, was 
prepared under Contract No. MDA 903-80-C-0196, entitled 
“An Analytical Model for Developing Objective Measures of 
Air Crew Proficiency with Multivariate Time Sequenced Data." 
This second volume provides documentation of the basic 
computer programs used in the processor termed "“Measure- 
ment Analysis Processor" (MAP). 
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The computer programs were designed to make ex~- 
tensive use of subprograms so that common functions used in 
the various programs could also use common subroutines. 
Furthermore, the structured approach to programming was 
used - to the degree that the unstructured FORTRAN language 
would permit. However, in the initial programming efforts 
limited use of the "go to" statements was maintained; but 
it was necessary for program variation to insert some "go 
to" statements which of course complicated the program 
structure. We hope this additional complication will not lead 
to difficulties in using or modifying the programs. 
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This volume contains programs used to analyze the 
subject data, procedures for use of the programs, and oF 
example outputs. 


On the procedure page, the user inputs are examples 
from data files at PMA. 


The computer used is an LSI/11, utilizing floppy discs 
containing programs and data. The DK: referred to on the 
procedure page corresponds to the disc drive used by the 
machine operator. 


Ihe subscripts indicating footnotes on the procedure 
page do not appear on the screen as computer outputs. 
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© NAPE INPUT FILE 
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Program Listing INTHEX< 
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Program Listing INTHEX (Continued) 
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Program Listing INTHEX (Concluded) 


Subroutine INTHEX ; 
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FORTRAN IV 

eeei SUBROUTINE INTHESCTHTAB, INTURR, JHEX HEK2D 
esez DIMENSION IHTHE(257), IBITSC16) 
cans TFCINTUFR.ES. ~32768) GO 10 6G 
eens IBITSC1E)=4 

GBs IFCINTURLT.8) iBITSCI6)=4 
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BS 108 CONTINE 

19 IBITS(12)=@ 
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Table 1 


INTHEX< 


, Integer to Hexidecimal 


Purpose: To convert a data file to thexidecimal notation for 
better transmission from time-sharing facilities. 

Note: The procedure will change with different time- 

sharing companies. 


Program Listing HEXINT 


FORTRAN IV HeIF-1 = FRI. @1-RY-B1 68254: 07 PAE O81 


C PROGRAM “HEXINT FOR‘ 
C PURPOSES CONVERT HEXIDECIMAL DATA FILE READ FROM DIALCOM 
C BACK TO BINARY UNFORMATTED (PROGRAM INTHEX AT DIALCOM) 
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Program Listing HEXINT (Conctuded) 
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STOP “END GF FREAD HEXINT. FOR’ 
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Subroutine HEXINT 
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SUBROUTINE HESUHTSLaTAe. LHER, INTVALD 
DIMENSION IHifi< 15.5, HEX) IBITEC 16) 
LOGICAL#! Hic. cHEX 
L=17 
06 46 I=1.4 
DO 1@ Js! . fe 
TFCIHEACL | EGA) OO TG ae 
108 CONTINUE 
Fié 
28 CONTINLE 
JME=! 
DO 3G fKe2,5 
ta-! 
IBITS(L=1M ae | ALEK? 
38 CONTIN 
408 CONTINLE 
INTVAL=@ 
00 568 M=i.IS 
hat! 
TFCIBITS(M) ED. GTURLEINTUAL#Z¢ rit 
388 CONTINUE 
JRCIBITEC168 80.22 CHTURC=<-INTUALS 
RETURN 
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Table 2 


HEXINT.FOR 


Hexidecimal to Integer 


Purpose: To convert a hexidecimal data file to an 


integer data file. 


Computer Output User Input 


e@ Run HEXINT 


<CR> 
Pause--I/O data disc in unit 
one (DK:) and hit return. 
<CR> 
Name input file (assign to DK:) *DK:TTEST3.DAT 
nh. <CR> 
Name output file (assign to DK:) *DK: TTEST4.DAT 
<CR> 
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DINENSION JDATAC 256) 
PAUSE “SNSERT COMPUT DFTA DISC IN UNIT (HE COKE) PD HIT RETURY 
WRITEC7, 7820) 
CALL FESIGNCS.” 6-1. °NER “NED 
WRITE(7, 7068 
WEROCS,, 5660: FSRICT 
Imi=e 
Nee 

108 CONTINLE 
WITEC7. 71: 
PEAS. 5606 TREE 
TFCIFINGH.E: . ! 
PRISE SES) CHP GAT DISC IN WRIT ZERO CS) ARE HIT REEH 
WITEC?. Ci. 
CALL FES]. Th RO NOD 


4 “ Program Listing FMERGE 
| £ { FORTRON TU HBA FRI 1-181 0859048 PACE. Gt 

a ¢ PROGIN “FIERGE.FiR’ 

e rf C PURPEES HERE C-50 DFTA FILES BY STUDENT MUPRER 

3 C COMMMACT C38 
oe COON RESEARCH INSTITUTE 
oe € FORT RUCKER FIELD UNIT 

C 26-ALY-88 
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206 CONTINE 
: FEROC2.EM= 8 UN TOUR TDUPG. TIDUS. IDLE, IDG. ~~ 
‘ | JORTACD) , 121. Aa: 
ININ+! 


WRITEC7.79G5 i CAD DATA) DATA 2Z? 
JECCIDATACS FE. 2ST) OR. DATA ZZ? HE, SILT) GTS AG 
WRITECS) TOURS 0) TOUR IOUS. TOUS. 150th, 

i AD EL A 


IOUT=](01+! 
@ 10 28 
36 CONTIME 
CALL CLOGE.2 
TO 1: 
98 CONTIN 
CAL CLE. =: 
WRITEC?, f at i tH, 1a 
STOP “EME ut PROG: EMbkut FUR 
328 FORMAT CIT: 
3068 FORMAT C7: 
7O1G FORMAT: HET PUT FILE AGS UG ass os 
7011 FORMATS’ ENTE G ZERO) TO READ ARGHER FILE 
i “CR ESE § COMED TU CLORE GURL FILE 
7O26 FORMAT C’ HARE LEST FILE CRBSIG! TS bh 
7088 FORMAT THN kets 127 
1 YUE tess 7,97) 
7060 FORMATC’ ENTER 202555 HUMBER CD; 
718 oo Wan DoePEem Ee 1.99 10) 


om oes 


$39 88 FASARASSERRAS RAGS SSSSSRASSR RRA SARE 


Copy available to DTIC does pot 
1 Heemit fully legitie reproduction, 


seme ss 


Table 3 


FMERGE.FOR 
rE 


Merge Data Files by Student Number 


The program FMERGE is used when a subject's data 
is located in two different areas (i.e., 2 different 
discs). FMERGE will bring both sets of data to- 
gether to form a single subject file. 


____Computer Output _ 
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@Run FMERGE 


Seong 
Py ; 


Pause ~ Insert output data disc 
in unit one (DK:) and hit €R 
RETURN 


Goend 
® ¢ 


Name Output File *DK:TTESTS. DAT 


CR 
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Enter Student Number ay! 
(eR 
Enter 0 (zero) to read another 
file or enter 1 (one) to close 
output file 
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Program Listing DFLIST 


W HOM 1 FRI GS-1 Gls 3es47 a 


C PROGRAM “DFLIST F(A’ 

C LIST THE C-5@ DATH FILE 
DIMENSION X( 24. 10) 
INTEGER X 
HELOK=6 
PRUSE “UNSERT (ATA 0:5C 1H DK: AND HIT RETURY 
WRITEC7, 7088: 
CALL AESIGHT?. 6 0-1. ROW NOLS 
WRITEC7, 7A: 
REROC7,. ONL: UKTI 

108 CONTINLE 

CALL JNANICN AHOOCK MA, 23 
IFCHREK. EG. 2 OT 999 
DO 20@ }=!.24 
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7006 FORMAT RAE G4 BLE: 
7OIG FORMAT in MUNRER POR LOSE ons 
7011 FORMAT CI, 
7026 FORMAT dl iL fish? 
END 
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Subroutine INPUT 


WV HOlfe] FRI CIA ~S1 61316554 PRE 1 


SUBROUTINE INPUT CH RELOCK. MARKLND 
DIMENSION CALAN SG? 29,10) 

1P=@ 

IPL=@ 

00 16 Jet.i6 

DO 1@ }=1..4 


REAXHEMBeSo) Reed > TMD NAR. IDUNS, IDR, TDS, JOU. 
! (OATH, 321, 2a) 

WBLOCK=NEL id t-+! 
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C MAIN LOGP 
00 1a i}, fi. 
K2K+2 
C STUDENT NUFEEF 
XU, D= Oat: 
C PEKIOO 
KA? 
XZ Dam Oh i te 
C SE@EN NUMEEI 
K>+10 
OG & I=1.: 
K=K+! 
HTKAOKGAT St 
XG. DSR ah, 
DATACE UATE or 
28 CONTE 
K-12 
3B LPSHIA: | SHI? 
WALT bok Th 
IP=IP+Ae. 
@ 10 7% 
@ IPLAPS 
X(4, LFS 
CSTIQK XK (ROU SHES ALING 15 INCHES TIE: 106 ites 
0 % 1-4, 
K=K+] : 
VARLORT DATA! cl SO STGS UALR 188 
XA D=nT’ 
3B CONTINE 
C STICK 9 (PITH GAS. OGG 1S INCHES TIES lee Miz. 
OO 6 $=9.1. 
K=-K+) 
VPLCAT DAC OL SQA SIGS UAT e168 
XC J=INT Sepy available 
@ CATINE pernit fully Ie 
7 Oe a 
Keke 
C SIDE TRY “Tis... 
NIK 
ORTACK e.. 


pesteneny en; cite ge 
ee | 
‘ i 


{ 
tI 
C 
c 
E 
C 
£ 


14 
15 
016 
17 
18 
19 
ze 
2) 
822 
23 
0824 
25 
0826 
027 
onze 
oa38 
Wi 
0882 
Ss 
04 
oS 
ea 
7 
ese 
wS9 
As 
i 
2 
en 
gal 
05 
we 
047 


to DTIC does 
gibl> pieducton 


omen I ce EE os 


ee ae FRE AE I SR eos ae : 


~~ ~~ 


agagna 88 § §§ 55 95 §8 $2 2 4% FE FROEGH BAAIAA 89S 


Ae wee ge ple nest 
net amen noe) 
. . * 4 @ ’ 


ees 


Subroutine INPUT (Continued) 


Vv WOlA-] FR: US-HAY-G1 GLE PRE @BZ 


X26. JAK 26, +1 OTK 
C PEDAL POSTION SCALING 15 INCHES TIMES 160 M3." 
YaFLOAT DAT ACES 40, 00951 S41 SIGH(S, DATAKK)9/16,6 
KCL JDSINTC 
76 CONTIME 
C PETCH INDICATED MAGES DEGREES, SCALED THES 106 
00 & I=15. 16 
KxK+i 
DATACK=DATACKO ADO SDATACKD 2) 
YaFLOAT CDATACE) 346, S49C5¢1 SIGHS, DATAKKS 18, 
XCLDSINTY: 
68 CONTINE 
C ROLL INDICATED Mii=3G CEGREES. SCALED TIMES 166 
06 9 I=17. {6 
KK+! 
DATACK Ceti} 4 SORTA, 2) 
niet Yoo 4902415 ONS, DATAKK) 971.8 
XCLJSIH7 


K=X+1 
X19, Doha. 
CTORQE 6 TO sit 
K=K+! 
X28. D0 % 
C AIRSPEED MAC=i4:..° 
XZ, Daeeh 
C HERDING 6 TG [4G {t- es 
K=K+} 
X22, FLO * 
¢ VERTICHL VELCC IT 
Kx} 
(23, D=0K 
CW BTC 14.5 (te. 
K=K+) 
X24 ELA Ce bots etedoe) a Slat dete 
C COURSE DEWITT IC. 
K=K+} 


MAED THES iit 
56 
TORLED TIMES ica: 


rie Aes ROT eee ees OPT 


Sami TIMES 16 


XB, stuns 
C SIDE ThEK NLbtti. 

KX-2) 

X26, D=f0 
C GENERATOR FEC‘ i: 

X29, D=OHn * 


See ati sadly Tine Idee 


e 


tind ifamki Ae. 


ue 
‘. 
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Subroutine INPUT (Concluded) 


[ FORTRAN IU HOlfr) FRI GitH-Gl O11 


spa es 


Dre ea 


wes ETI 

I 9069 681 CONTIN 
28 RKO 
9851 RETURN 

ri 2 002~COt«C«éD 


FBR R ES Ue 


SO 


i 
i 
4 
L 


me aS 


~ om 


' 


; eo5 s 


| rey iss! Owe! 
e ’ e ‘ ; 


o~ 


I 
i 
I 
'§ 


Awe ‘ e ee | 4 e ‘ ‘ 


Table 4 


DFLIST.FOR 


SSS eee 


List of Data File 


Purpose: Reads the data file and presents the information 
on the CRT. 


Computer Output User Input 


e Run DFLIST 
<CR> 


Pause -- Insert data disc in 
DK: and hit return 


Name input file *DK:F 20935. DAT 
<CR> 


Enter unit number for listing 
(6,7): ‘ 


line printer 


teletype 


4 gh ic car ead ae a 
SL watt a Parkin hos 


7, * 
oes 


Example Output - OFLIST = 


~~ Fs om 


BIS AAS PAS AS FAS OSS ATS Subject No. 
‘ ‘ Session No. 


“MB -N3 -73 -AR IRIS. 16-16-57, -18R_~—- Stick &X Position 


¢ 


Pe a ee | 
¥ t 


“MM BW MW WwW MW M BM MW MW Stick 4 Position 


ner a AR an = Se 
. 
[ere | 
* ; 
> eb oh tp OD ee Oe EP ee SP ee 


] Vv oT OURO /—COKCOC-BH™s*C«=%EE”—«~PP@daal: PPcsitiiorn 
eis “% MN -% -% -% % % *% -% -® 
we GF MM MA Mm ss 8 WR We 2B Pitch Angle 
\ oe MG MH MRM 6 MRM RW 4H Wir 
ae “KB KB IN? -164— -184 -1R6 - 201-210-2108 --245 ~Roll Angle 
“7B NGS -1514 -1847 -1859 -1912 -2108 -71% -Aie -2%. . 
: “- WG NSSL 97TH 1T2?—sN79s1A OEE Altitude 
vie 5 S S BS BF BS BF BHR BS SB Torque 


38 MH 6 Mh hh lM? MIs Airspeed 

1KA2 IRA? 14916 N4NS 7A INEM 2 IT fa7is | Me Heading 

$012 1178 14 Re SH ? Vertical Velocity 
G6 @ “BW NN +6 46 eR -R -%H H OV aw Angle 


~ OF 


BSBNSRORBRNSLRBSSSRATASHSveveuaun~ 


+4 


‘ 
‘ 
® a 
ai! 
1s 
1 oy - 
fw 
{ wit 
ye # 
* 


4 
Seen 
° ‘ 


i ot 
available to DTIC does n 
pose fully legible reproduction 


Saeed 
eaeey 
é e 
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eee 


id 
| F 

‘¢ 

tL 

| | 
: [ 
| 
be 
: ais 


Program Listing STRPLT 


| W HOlARS «—-FR] GR-+HAY-21 GR e589 PRE GS! 


C FAIRPISE 15 1G PLOT (-5@ URRIABLES ON LINE PRINTER 
COATES LY 2. 1968 
C ARON RESEARCH INSTITUTE. FORT RUCKER 
i LOGICAL+! LHEAD. LPCHAR. LYTIT. LYSCAL, LXSCAL, 
i LYCINE, LAL INE -LUARC, TESCAP. LTY. LYVAL 
WBZ DIMENSION LPCHAR( 10). IBORTAC29, 16) ,LHEAD(68). 
. J LYSCALC1@1>. LYLINEC1@3), 
i LUPRC (16). FRO 16) LITITCES), 
i LYVAL(104), INDEX 16), XMINC16), SHRKCI6), 
i RANE 16) 
DATA LXLINE/124/.L LIKE -1a3eq5-, 


i C PROGR: “STRELT FOR: 


: «, ’ ¥ 
g ; 


1 LYSCAL/126.$052. 126. $432, 126, 9032, 124,932, 
1 124.9432. 124, 9932, 124,9832, 124.9432, 
! 124.9932. 124, 9432, 1247, 
. 1 LYWRL/°@ 9920.1. 0" 8832, 727, 70" 8932, 
! 7B BHTD. 4" 787 BHR, 5, R8032, 
1 ! “BW BASZ 77 OY £8032, “RO 8432, 
4 | ‘ 1 6G: Te Be $70". 7 “ft, 
$UMUIT ASC RRS PE 
a ! bol OR OWT B82 
ue Soe | 1 LURRC 02 8.5.11. 114,97, 
5 ee i O12 W122 AZLHA27, 
; 1 ESCRRS 27. LTV ISP TU feed, 
a I 10 IMME 5.4.13 $8.17. 19,20.29,22,23.24. 540 
is id PRISE “PEADY LIKE PRISER AMD HIT CORTIMLE’ 
é 5 PREE UHGEPT CaTa C120 1K UNIT (HE COVE APD HIT PET 
& 16%. MPYTEC? 7a 
®. Er, (AL ASSIGEED 0 =f HED 
(8 UPITEC7 PG: 
19 PEAS SEI? LHEAR ST TET. 68) 
4 WITES 7.7 
% PEACE SE MESS Deh AD 
WPILTES? OTS 


EROS SEE) FAMICT 
WPITE 7 BC. 
PEAS SAR TORT? 
WHITES? A 
REPOS SGP: THEGTH 
CHL SCALECIROATA iH. SRY MAP. PPE THES ISCICTS 
IF CNPEA ECG) Gf) Tih Oats 
WPITES HTT Eiie> 
WPITECIUNIT cea) ‘LHEAY 1). 151.68) 
UPLTECIURNT CA 
OO 12 =f! 
TECIUMUMCD) MELD) Of To 1 
WRITES TUNIT 6046" 1 INDES DU NINCD RECD RAKED 

10 CONT IIE 
WRITE TUNIT 6080 TD LURPOCD) Jet. $1) 
HLOK=@ 
ILOC=? 
ISNT =O 

"108 CONTIMLE 

CALL SMPUTC IEC Cre Rey! 29 


19 


PCR T RE Tt! HOIA-l = FRI GO+RY-O1 68:65:69 


AL AIL GREAT So Tene nty tt 


5 
4 


Program Listing STRPLT (Continued) 


IFCHPRK.EG.@) GO 10 88 
IFCHPRK.E@.-2) @ 10 818 
IFCNBLOCK.NE.1) GO 10 128 
WRITECIUNIT 6418) 

00 115 I=1.16 

TFCIBDATAC!»1).HE, 1SBICT) GD TO 115 
00 11@ JRL. 13 

JSIBRINDER(S) 


WPITECIUNIT 6408) J. ISB, IBDATACISIB, 1) 


HG GUATINUE 
GO 116 

115 COTIME 

Me CONTINUE 
bRITECIUNIT 6851) 
IFCIUNIT.EQ.7) GO TO 128 
UPITECG, 6658) IESCAP.LTY 
WRITE(6.6878) (LHERDC1), 1=1-68) 
UEITEC6 6408) CLYTITCI), 1=1, 108) 
WITES6.6298> (LYVALCI). T=1. 164) 
WPITES6.6108) <(LYSCALC]), 121,161) 
WRITEC6.61100 CLYLIMECT), 121,163) 

{26 COAT IBLE 
10 308 I=1,16 
JFCIROATAC1).HE. ISBICT) GO TO 388 
ISPRNT=1SPRHT +4 
TALOC=1XLOE+1 
OG 2728 J=1,108 
LPCHARC J)=5, 

72 COTINLE 
1X6 268 Kel. $1 
TFCIUMUM(K) NE. 1) GD TO 268 
ISUB=THDEX(K) 
TUALUE=IBDATACISUB. 1) 
JF(JSUB.E0,22) IVALUE=1ABSCIUALUE) 
PURLUE=ABS (FLOAT ( TURLUE)-YHIN(K)¢!.@) 
YLOC=RUALLE/ RANE (K)4108.66.5 
JFCIUHIT.E@.7) YLOC=YLOL/2.@ 
JYLOC=IFIXCYLOL) 
JFCIYLOC.GT.108) IYLOC=108 


JF(CIUNIT. EG. 7). FAD. CIYLOL.GT.58)) ILO 


IFCIYLOC.LT.1) IMLOC=I 
IFCIXLOC.LT. IREGIN) @ TO 738 
IFCIUIT.E0.7) GB 10 238 


WRITE(7, 7308) NBLOOX, 1K, IBDATAK4, 1), TUALUE, TXLOP. IYLOZ 


728 COTTE 


TF CLPCHARCIYLOC) £6.32) LPCHARCIYLOC)=LURRC CK) 


IFCIAIB.HE. 22) GO 10 268 


1 LPCHRRCIYLOC)=43 


IF(CSHOME.E.6) AED. CIBDATAC ISB, 1).L 7.8) 


1 LPCHARCIYLOD)=45 
THSAVE= 1BDATAC ISB, 1) 
268 COMTIME 
LESCAL=32 


4 
oe [ 
t be 
| ‘ : 
i Program Listing STRPLT (Continued) 
FOR TM 1 HRIAEL «FRI @S+IRN-81 @ase5:09 PAE 083 
3. aes IFCISPRAT.EG, 18) LESCAL=45 
tees IFCISPRAT.EG.1@) ISPRNT=2 
meg IFCIRLOC.LT. IBEGIN) GO 70 380 
i ett! IFCIURGT.EG.7) GS 10 278 
» ais WRITECE.6128) IBOATAC3, 1), IROATAC2, 1). IBDATAC4. 1D. 
1 LYSCAL-LYLINE, CLPCHARCLIKL). 1JYL=1. 198). 
a 1 LOLINE, LXSOAL, IBDATAG4. 1). TSLOE 
| a. O16 @ 10 3 
M15 | 27 CUTIE . 
1 -  Ots WRITE(7,7128) IBORTACS. 1D. JBDATACZ, 1). IRDATACG. 1) ISLO, 
{ = 1 LASCAL,LXLINE, CLPCHARCIJEL). WML=1,58), 
. he 1 LSLINE, LXSAAL 
W117 32 COMTINE 
/ ANE i To 108 
~- G14 | 0 CONTIME 
i em WITEC6. 710) 
| ates 6 10 8m 
5. WZe | BIG COMTIME 
we IFCIURNT.EC.7) Gt TC 828 
WA WRITEC6-61IG) CULIHECD)J=$.403) 
aie WRITEC6. 6108) (LISCALCD). 121.180 
a? WRITEC6 6890) CLYVAL CL). J=1. 164) 
aly WRITEC6 6081) CLUTITCD, Tst. 108) 
alr WPITEC6-C270) CLHEANCT), 121.68) 
3G 82 COMTINLE 
ez! CAL CLOSEC?' 
«137 STOP “EMD CE PROGRAM STRPLT.FOR’ 


HSS | AIG FORT CORAL) 
O124 SEG FORMATCHI17) 
M125 SEES FORMATCI?) 
G136 S35 FORNATCI7) 
GSS" S63? FORMAT CI?) 
@32%  7GGS FORMAT(’ NAME THE INPUT FILE CASSIG! 10 [k:)°/) 
@129 POLS FORMATC’ ENTER TITLE FOR PLOT <6@A! MAX)”) 
C146 PETG FORMATC’ ENTER FLAG FOR THE ELEVEN URPIABLES (1117)") 
G1d1 7855 FORMATC’ ENTER STUDENT NUMBEP <I>”) 
W152 73S FORMATC’ ENTER GUTPUT UNIT HUMBER <LP=%.1T:=7)") 
C142 7037 FORMATC’ ENTER BEGINING LOCATICH OF PLOT ¢1)’) 
G44 710K: FORMAT (’ READ ERROR IH INPUT FILE") 
(145 72S FOPMATCIX, 717) 
C145 6618 FORMATC’ 7.5%, “C56 STRIP PLOTTER ROLTINE’- 
1 6X, “ARO RESEARCH INSTITUTE, FORT RUCYEP FIELD UIT’” 

Pa 1 6X, PULTIURRIATE ANALYSIS OF TIME SEQUEMDED DATA’ / 

: 1 6X, "UH-I (HEY) STUDENT PILOT SIMULATOR SESSIONS’”” *) 
. O14" 0626 FOPMATC 7.5%, °TITLE: “.60A1) 


oe 


enemy 
. ' 


MS GEOG FORA’ “7° “5%, “URPIRBLE SORLESS“7° “7 Copy availabls to DTIC does not 
: 17 7.5%, /INDEX’, 3%, “SB’, permit fully legible reproduction 
i om! MINTY 25, “HERIINI 6%, “RANE” 
f “ 1° 7S, 


! L Lhe a Re’ 
CSE FORT CON IS, 24,14, 2%,F7.8,28,F7,8: 20-98) 
CEE FORAT CIR Als 1°) 
C151 ASG RIGT IS ff 

1°50, URPIRBLE REVS “77 “7 


eo, 
4 e 
~— © 
=o 
cn te 
43 


21 


nm 


; 
} 5 
n oek 


‘ 
: 
H 


i 


FC 


3 bad 


HEI ARI 


€X, 13. 20 AL SE, 


OH 1S. 28. AL SK 

Sy QUAL SE, 
O85. 28. AL Se 
68,75, 2 AL SY 
OX. 2K AL 
6815. 2K AL 
NTS Qh AL 


SR, “TNDEX 28, “CHAR” 2K, “DESCRIPTION / 
hy Ky << 


Program Listing STRPLT (Concluded) 


FP] GS-81 83:65:69 


“STICK X CROLL INPUT) “7 
OY. 15. 2K ALS” 
6M. 352K, A1 SY,“ 


STICK VY (PITCH INPUT)’ 
PEDAL PRSITOH’/ 


“PIICH’? 

POLL’: 

ALTITUDE’ 

a TOMES, 

ts AIPSPEEDY 7 

*s HERDING’ 7 

‘ VEPTIOAL VELOCITY: - 


tet fy 


PRE C54 


H 
i 
! 
1 eM! 
1 6K] 
1] 
{ eI] 
1&1 
16h] 
Le 
1 6%] 
! 
J 


CHATS. ZALES. WR Eo) 
BEST FORST YA e. HOE 1: ALL SORLES PELATIUE EXCEPT HEADIN: - 
1 6S. NOTE 2: URPIABLES SAMPLED OPE ELEY Th SECS’ » 
1 6. HOTE 3: SVPBOL "+" INDICATES HEADING HG’. 
1° CHASED FROM NEGATIVE TO POSITIVE’: 
1 6X. “HOTE 43 SVMBOL "=" THDICATES HEADING HAG’. 
1° CHRAGED FROM: POSITIVE TO HEGATIVE’<” 7°) 
CES FCOTAT I ORL IR IX) 
BLS FORMAT IGN, 100817 1%) 
CESS FOOHAT CIR. 108AL 18) 
BOG FORMAT CITE, 10481) 
BGC FOPMAT (17H, IGIAL) 
BLIGE FORMAT (16x, 1634%) 
S128 FOPMAT (4X, 12, 15.14, 1%, 10401. 14.76) 
7128 FOPMATCAX, 13.1%, 15, 1X, 15.1%, 15,4, 4A) 
6416 FORMAT COX, “IHTIAL UALLESS “71X/ 
3 6X, “INDEX 2. “SUB” 4X, “URLUE’’ 
1 6X7 AK eer’ *) 
C462 FORMAT (6X, 15, 2K, 13.2%, 17) 
END 


S55} 
fe’ 


(S835 
As4 
er 
an 
ke 
sx 
> (NG 
mi: 
ars 
(tS 
ete 
Cat? 
G19 
7 
3 


3s, 


“FOOT IM! 


Subroutine SCALE 


HBIARL «= FRI: GO+RY-81 08°045:57 PRE C2! 
SUBROUTINE SCALECIBOATA. XMIt. SHAK. NAPE PANE. IHDEX. 158 1CT) 


DIMENSION IBDATAC 29.18) OMINCI6) SHR fe). 


! BIRCIG) PRR CIG) RAEIG) THER 16) 
HELOLK=B 
06 16 I=f,f6 
HNC) =32767 
MAK CT =-32767 
1 COMTIWE 


1 COMTI 
CALL INPUT< IBDATA. HELOMY’ MARY. 2) 
TERRA E10) Gt TU 46% 
(Ct 2H I=1,1¢ 
JFCIROATAC! 1) HE. ISBICT) GA TO 288 
MW $66 Va. 41 
DBATDEXKK) 
IFCISUB.EC. ZZ) IBDATACITR. D=1ARS<¢ IBDATACISB. 1) 
IFCIBOATACISLB. DLT. MINK) MINCK=IBDATACISUB. 1) 
IFCIBORTACISUB: 1).GT.MRAKKK)) MREKKOSIBDATACISNB. 1D 
156 COOTIE 
208 COAT INGE 
(i TO 168 
4 CONTINUE 
WITEX7, 7688) 
00 56 J=111 
WRITEC?,7028) L-MINCL) MARCI) 
Sis CETIME 
IX! 668 J=1,16 
SHIN T=FLOBT CHING I) 
SHREK D=FLORT MRED) 
TFCCMINCL).LT. >. FRD. CMARKET ).LT.@)) 
1 PREC TD =ABSCOMINGT) FRG CHAR 141.8 
TFCCHIFKT) LT.@). AOD, (MARCI) GE.@) 
1 RANGEC D=RBSCANINCD) 40K 41. 
JFC(CMING]). E68). AND. CMPD) LT.@)) 
1 PRE CD =ARS CORR 1) 4X87HC 1941.68 
IFCCMINCT).GE.6) AAD, CHARC I). CE.) 
1 RANGEC 1 =XMARCD)-XMINC 141.8 
6698: CONTINLE 
MMIC 9)=8.68 
KS) =19088.€ 
PATE (9)=1908! .8 
PEUITD 2 
PETLPH 
7G FORTAT(’ URRIABLE MINIMUM FRD MAKIMLM URLUES’> 
7028 FORMAT (3X, 15.2%, 17.2%, 17) 
3, 


Subroutine INPUT 


Foote y' HOIA-] = FRI @BHARY-61 08:05:53 PRE G2! 


eet SUBROUTINE. INPUT(K, HBLOOK, NRK. ND 
ere INTEGER DATA, X 
(iss DIVERSION DATAC258) -X(29, 18) 
ae SS 1P=@ 
(x5 1128 
(Hs 00 18 JEL. 18 
eee 00 18 121,29 
Cees KL I=@ 
(0413 CONTIME 
cst tre =8 
nit PEROCH ENDS6it -ERR=O@1) TOUR. TOUPZ. TURE. LOL. TOURS CUP. 
‘ ! (DATACT >. T=. 258) 
eat: HELOCKIBLOCE+1 
e153 Kee 
{_ ¢ NATH LOOP 
eats 00 108 Jet. 16 
eats ka4e2 
¢ STUDENT HUPBEP 
aie XL DRT 
¢ PERIOD 
"ay? rKe2 
an G2 X(2, J =HOD(DATACE) .256) 
( SEGENT HUBER 
ents KxH1@ 
wan 00 28 11,2 
at KK! 
ee NIV-=NOCDATACK) 2) 
ee 13,.D=K(3, DAF 
24 DATACK=DATACH) 4ITF 
1 2 CONTIN 
(THE 
26 KL2 : 
(OF 3 LPSHATACK “2564 IP 
enc IFCIPL.LT.LPS) G0 10 48 
ws IP=1P+256 
wrt OIG W 
OZ 4 IPL=LPS 
53 4, J)SLPS 
( STICK CROLL INPUT) SORLING 1S INCHES TIMES 108 MEDS. 24" 
ea4 158 15.8 
acs Kes 
PET YSFLOATCDATACK))#8..01904441SIGHCS, DATACH)) 16.0 
en7 HCL, JSINTCY) 
(2% «SE CONTIN 
C STICK Y CPITOH INPUT) SCALING IS INDHES TIPES 108 MAE=6.333" 
(4 00 6 19,12 
(40: Keke! 
(41 Y=FLOAT CDATACK) 40. 0193234 SIGS. DATACH) 18.0 aa 
47 XL, JDSINTC) Copy available to DTIC does i 
(4. eo CONTIN pecmnit fully legible repro 
4s 00 78 1213.14 
Os Keke 
(SIDE TASK STRPT FLAG 
(646, HIV=ANOCDATACH .2) 


gic 


x 
% 
u 
$ 
a 
te, 
is 
z 
& 
& 
&. 


I 
I 
i 
I 
i 
[ 


fay Gren. 5 
. . 


Yaa ae oe ee 


_ FORTRAN 1U WOIR-! = FRI CB+HAY-6! 08:08:55 PRE O82 
47 DATAKKOSDATACK) TV. 
es X(2B, 38X26, J) +1 CNTY. 
C PEDAL POSTION SCALING IS INCHES TIMES 108 M=3.25° 
49 Y=FLORT COATACK) 198. 6899194 1SIGK(5. OATAKK) 16.6 
iwi] KCL DDSI CY) 
5! 76 CONTIME 
C PITCH INDICATED MAK=66 DEGREES. SOALED TIMES 108 
Cz 00 88 I=15.16 
A: KX+1 
4 DATACK=DRTACK)-HODXDATAKK) » 2) 
55 YFLOAT CDATACK) 90. S495 1 SIGHS. DATACH)) “12.6 
” 5 XCL SINT CY) 
005? «= CONTINE 
C ROLL INDICATED MAN=98 DEGREES. SDALED TIMES 13 
0058 00 9% I=17,18 
9 KX¢+l 
aa) DATACK=DRTACK) 0D (DATAKK), 2) 
O41 Y=FLOAT CDATACK) #6. 54932415 1Gi(S, DATAKK) 718.6 
(62 XCLJSINTCY) 
53 «9-9 CONTINRE 
C ALTITWE 
0064 =K+1 
0R55 X19, D=DATACK) 
C TORE = TO 108 PSI 
RES K=Kel 
0067 X28, J)=DATACK) 256 
C ASPSPEED PIR=146.5 KNOTS, SCALED TIMES 168 
ORES X21, J)=6OCDATACK) . 256) 
C HEADING @ TO 368 DEGREES. SOPLED TIMES 1@ 
(BES KK+l 
7 X22, J) =FLORT (DATACK) 98, 5495.24151GK5,DATACK)) “10.8 
C VERTICAL VELOCITY 
wri KX+1 
OFZ X23. D=DATACK) 
C YA @ TC 14.9 DEGREES, SCALED TIMES 168 
ea BOS! yap LontcoaTAack) 296. 8¢54R+1SIGICS.LATAC)/108.6 
care © COR AGUA 
6 X25, J) =FLORT CDATACK.) 98. 3493.24 1S 1G5. DATACK) 16.8 
C SIDE TASK NUMBER 
Or? K-21 
e738 X(26, T)=K(26, 34D CDATAKK) , 256) 
C GENERATOR RECYOLE TIPE 
C79 X29, JSDATACK)/256 
C STATION SELECT TIVE 
L028) 2K+22 
ee! X(27, J2DATACK) 256 
C NIE RESPORSE TIME 
SZ X28, J) ANOCDATAKK) » 256) 
= «= 108 CONTINE 
(924 Fe =18 
(BS PETUPAI 
Oe | EAL CONTIILE 


Subroutine INPUT (Continued) 


Subroutine INPUT (Concluded) 


PRE 083 


FRI GO-HAY-6! 68:65:53 


HatAR! 


} 


Table 5 


STRPLT.FOR 
SSE 


| 
Strip Plot 


Purpose: To plot flight variables on the line printer. 


Computer Output User Input 


ae ai AiG! 
Pause -- Ready line printer 
and hit continue <CR> 
uni e i i 


Name the input file (assig *DK: F17932. DAT 
to DK:) CR 


Enter title for plot (goal Max) Subject 17932 


Enter flag for the eleven 0,0,0,1,1,1,0,1,1 ,0,0" 


variables (1117) 
<cR> 


Enter student ‘number (1) 


= Unwanted variables 
Desired variables 


°o 
I 


Table 5 


STRPLT.FOR (Concluded) 


| ? Enter Pa ieee unit awer 
(LP: = 7/8 


Enter beginning location of 


Pree] 


‘ 


ay hy FY be fF: 


omng 
es: 8 


—_————— PR 


oe tell Oh de BN em ee 


ay i YY 


igthegeey meg orsst aS ORE * 
‘ e . 


ERR Ee 


(-56@ STRIP PLOTTER ROUTINE 

OO’ RESEARCH INSTITUTE, FORT RUCKER FIELD UNIT 
RUTIURRIATE APLYSIS CE TIME SEQUENCED DATA 
tHET CHEV) STUDENT PILOT SIMULATOR SESSICHE 


THLE? ABIECT WSS 


VARI HELE SCALES: 
OES SUB OMNI MRP PANE 
4 §5 -35. 166. 1971. 
% 47 -1861, S21. 4813, 
6 49 Wi, | 14%. 506. 
€ “a 768, 97%, ai. 
“4 72 ¢, eae, 19001. 
WHRIHBLE KES: 
IES «CHR DESCRIPTION 
1 STICK ¥ CROLL INPUT) 
7 P STICK ¥ CPLICH INPUT) 
an PEDAL POSI TON 
4 Ff PITCH 
Sor ROLL 
6 a /ALTITWE 
rot TEE 
ay AIRSPEED 
Sh HERDING 
1a 7 VERTICAL VELOCITY 
fo wu 
INTIAL. VALLES: 
Wwe: eh = UALLE 
a) 112 
v 4 7% 
2 -38 
4 98 vi 4 
% 47 aS 
& 19 | 
Mi B 


“$e 

4m -t421 

mM B 20 

1 om 10 
MHTE 34 Co ACES SARC de Boe POD Gone 


HITE 22 SYMECL "4" INDICATES HERDING HAS CHANGED FROM NEGATIVE TO POSITIVE 
_ ME 42 SHPROL °-° TNDICHTES HERDING HAG CHANGED FROM POSITIVE TO NEGATIVE 


Example Output STRPLT 


Example Output STRPLT (Concluded) 
Y ARIS PERCENT 
A 


% .. 


! 2 
b Os oe 
i = 
F: aa 
J 
Fs 
: Wa * & as Posesce ie acne ees 
a I a va ar a ea ead & -|-« 


BEET HLS SLES RASRSE SEEKS OSE SLRS ' 


xrceeeee Tee OHH OSHS FES SRE RE SESE 


SSE ESF ESSE SESS SSS SSL SSSSEW SSE EG 


— on ow 2 S88 8 we SE OT 


eS ee Bi Sara care anal 
be ce OY Se 


ae 


= ee be FE! 


—— 0 lee pe 


6 


Program Listing PHAN 3 
WIA-3 TH G7-V-81 8:52: 19 PRE @8t 


PROGR: PHRIG. FOR 

PURPOSE: RETURNS CELL URLUES FOR C-5@ FLIGHT DATA 
CONTRACT: (58 

DATES 21-FEB-81 


COMIN URRELE (3088) » ERROR (3888) 
DIMENSION ICELLS(25), 1 THEAX(25,25) -KSTART (108), 
IKSTOPC 108) ERIK 4, 108) 
JNTEGER VRBBLE. SUBICT, HUPURR. NUPL SUBI  SUB2, SUBS, DEV, CELL CALLIN 
REAL 1,J,CHKUAL,ERIERN 
DAH TOTERR. TOTSOR, ERRSUR, ERMERN®. 6,8. 8,6.0,6.0/ 
DATH RATERK, TOTKTE RTESOR, ERRSOR/@.0,0.6,6.8,8.0/ 
DATH SQB1,SUB2, SBS, NUMJAR 8, 8, 4,0/ 
DAH ICELLS/ Awe’, TTRRTR 62500 


CIMT INE 
PRAISE’ HIT RETURN & NAME CUTPUT FILE’ 
CALL RESIGN,’ 7.71, °NEW’, “NO“.1) 
PHISE’ INSERT DATA DISC IN DRIVE ONE AND HIT RETURN’ 
PAUSE’ HIT RETURN AND NAME THE IST INPUT FILE’ 
CHLL AESIGNC?, ‘7-1, °RD0", “NO 1) 
WRITE(?7, 7822) 
FADS, S77 BIT 
WRITEC7, 7488) 
FERDG, SARA CHUA. 
WRITEC7 , 7872) 
READS. 5073 ERAEAN, ERRDEV, RTMERN, RTEDEV 
WHITEC7, 7010) 
FFADCS, SRA) NUM 
JFCMUMLFQ. 73) SORLE=10.8 
JFCNUPLER. 19) SOALE=1.4 
JECHURLEQ. 22) SCALE=108.0 
JFCMUM. FQ. 23) SCRLE=1.8 
WRI TFC7, 728) 
READS, 5628) NEV 
00 71 Ti=t,1@@ - 
EAD, 5821 KSTART (31) KSTOPC 31), CARER J3,33),J3=1,4) 
JECKSTARTCIT).LE.OG0 10 22 
CONTIN 
ILIST=3 3-1 
(ALL (L08E(2> 
PAUSEHIT RETURN & WE 20 INPUT FILE’ 
CALL PESIGK2, ’ *,71,°O0’, WO’, 4) 
DO 2 Viz, ILIST 
TF CHUN. ER, 19) ERNEROSRPERIK | 13) 
JF CRRLER. 21) ERREROMAPERC2, 33) 
JF ORM. EQ. 22) ERERMIERKS, 11) 
TFCRSLER. 73) ERERARERIK4, 11) 
TSHIRTS WRT 3) 
ISTOPSSTOPCII) 
CALL (EDI COAL ILORAR, ISTART, ISTEP, SBICT) 


ERR=CFLON CVPIBLE (K) )/SORLE)-CRUAL 31 


Program Listing PHAN 3 (Continued) 


WIRE] «= TH @7-N-81 68252: 19 PRE 082 


ERROR(K)=ERR 
CTNTIME 
NURRI=MURR-1 


UR) TECEY, 708) SUBIC. NM, 
1 NUNVAR, SCALE, ISTART, ISTOP, ERNERN, ERRDEV, RIMERN, RTEDEV 
WRI TECDEV, 7851) 


PRERMERNE( 5. SHERRDEV) . 
E-ERPERN+ CERRDEWHO, 5) 
(=ERPERN- CERRDEWHG, 5) 
D=ERMERN- (1. SOERRDEV) 
F=RTHERNEC 1. SORTEDEV) 
F=RTMERN CRTEDEWH@, 5) 
GAR TMEAN-(RTEDEWN@, 5) 
HER TERN (1. SORTEDEV) 


cette 


WRI TECDEV, 7855) 
WRITE CDEN, 7878) 
WR TECDEU, (61) 


LINES] 
DC 4@ HEI WURRI 
J=RRORH) 
J=ERROR CHL )-ERROR(N) 
(tu=@ 
IF(3.CE. END. 1.0E.A) CHL=5 
1F(3.GE.E. AND. 1.CE.B. AND. 1.LT.A) CELL =4 
IFC). GE. E. AOD. } 8.0. AND. J.LT.B) CAL=3 
IFCJ.GE.E. AOD. 1.0.0.0. 1.LT.C) CHL=2 
JF(J.GE.E.AOD. JLT.) CELL=I 
TF C1. GE. F.AND. JLT. EAD. 1.E.A) CELL=10 
FCI. GE.F. AND. JLT. E. AND. 1.08.8. AD. 1.!.1.R) CELL=9 
FAD. J.LT.E.AND. 1.GE.C. AND. 1.LT.B) CLL=8 
FAD. J.LT.E F000. 1.0. AD. TLT.C) CLL =? 
FD. JL T.E. AAD. 3.LT.D) ae 
F.AND. 1G.) CELL=1 
J 
} 
] 
] 


IFC. oF, 
IFCJ..F 
JLT. A 

1F().GE Sol. PD. Cmts CHLL=14 

FG. fe. EARL ae. o(E.C. AD. J.LT.B) CELL=33 
JF(J,.GE.G. AND. J.L1.F AD. 1.E.D AD. LTC) CELL =12 
TFC3.GE.G.FOD. J LT.F. AOD. TLIO) CL=1S 

1F(J. GEHL AND. J.LT.G.AND. 1.4E.A) CELL=28 

IF(J.CE HAND. J.LT.G. FO. 3 .CE.B. AD. 1.1.) CELL=19 
JFCJ.CE HAND. JLT.G. AND. EC. AOD. 1.LT.B) CL=18 
IFCS.CE. AAD. SLT... 1.0.0. 1.07.0) CLL=17 

IFT, CEH AD. .LT.G.AND. 1 LID) CELL=16 

TECLLL THAD.) E.4) (FLEAS 

JEU T.HAND. 1.0F.6AMD. LTA) CEL =24 

JEQLUT HAE. 1.0 AND. LTB) CELL=23 

IFCLLT AAD.) EDAD LTC) CELL =22 

JECT. HAND. LTD) CELL=?2 

SECCEL LEGA) PAUSE “ERROR IN CELL DETERMINATION’ 

WHITE COE. em) 41,0, CELL 

ICH LSCCEL L=TCELL ACCEL 43 32 


HOIAR 1 TH! 7-4-8! 68252319 PRE 68S 


IFCWLL1.2) CELLMI=CELL 
TECHLLT.2) GO TO 39 
JTWATX<CELLINA CELL) =] TMATX (CELL CELL 41 
CELLMS=CELL 
3 CONTE 
LINES=LINES+ 
JFCLINES.LT.48) GOTO 46 
L IHES=@ 
WRITE CORY, 7871) 
44 (INTIME 
REWIND 7 
BR Cie 
Jt « WRITECDEV, 7898) CICELLSCM).MEI, 25) 
URITECDEV, 7893) (CUTMATXCALND NEL, 25) .Me1,25) 
D0 98 HEL 
CH WRI TEC) 5824) CU TMATXCPL KD NED 2) 
CALL (LOSECH) 
CLL CLOSE (DEV) 
STCP 
¢: 
“ae FCQMATCF 3.7) 
“AIA FCQMATCIE) 
“@i} FORMAT(F 13.7) 
S01 FORMRT(217,4F8.@) 
SOR FORMRTCI6) 
1677 FORMATCI?) 
“408 FORAATC4FI5. 4) 
£474 FCRPRTC2514) 
7A FCRART(’ ENTER CHECK URLUE TO BE ANALYSED’ 
j ” AIRSPEED-ENTER 98.0877 
! ” PLTDTUDE-ENTER 2088, 08” / 
1 ” HEADING -ENTER 98.087/ 
{ ’ RATE OF CLIMB-ENTER @.08’/ 
7HI@ FORIRT(’ FOR THE URRIFBLE TO BE ANALYZED: “7 
j * AIRSPEED ~ ENTER 21’/ 
j ” RLTITUDE ~ ENTER 19’ 
! ” HERDING - ENTER 22’/ 
j * RATE CF CLIMB - ENTER 23”/) 
7008 FORMATC’ FOR CUTPUT GLP 7217 ’) 
7P01 FORMAT(’ ENTER THE START AND STOP POSITION IN FILE, ENTER ERMEAN’) 
7KO FORMAT’ ENTER SLRJECT NUMBER’) 
7A. FORMAT’ ENTER ERNEAN, ERRDEV, RIMERN, RTEDEU’) 
7168 FORMATO’ SIBJECT! “167 
i * URRIPBLES ’,12/ 
j ” MOBER CF SPRPLES 7,157 
§ * VPRIFOLE SORLE: *,FIG.SY 
j * START POSITIONS 7,15/ 
j * STOP POSTION: 157 
/ * WEAN CF ERRIR = “4F14.4/ 
1 * ERROR DEVIATION = 7,F14,4/ 
! ” WEAN OF ERROR RATE = ’,F14.4/ 
£4 ” DEVIATION OF ERROR RATE = ’,F14.4/) 
781 FORIRTC’ 7”) 
768 FORMATC’ CELL ULES 7 33 


Program Listing PHAN 3 (Continued) 


Program Listing PHAN 3 (Concluded) 


i HENeeAT WY Hele) TH GP-HRA-B) 0252219 


a 
a7 
We? 
as 


74 | AOS FORMATE I” 


1 Ves) 


Lad 
we 
PY 
Pal 
s 
* 


ree ea en 


if FRR’) 


PRE G04 


TTR FORMATO 7, INDEX’ SR. “ERROR 20, RATE CF ERROR ZK. CELL 7° 


TA) FORMATO“ 7, INDEX’ 5%, “FRROR’ «2K. RATE 0F ERROR’ 2X, “CELL 7) 


WO FORMATO! ISIN F9S: PLS. 2, 2X Td) 
TG FORTS CELL ALES * 


SPINS SEN 


pst 


OPTS Ie 


= ttt 8 


- ft tt 8 


“SE SOM TS. SRY 


, $+ 4 tS 


J 

J 

J 

J 

J 

$7 NIN TS. RATE OF ERROR. 
1 : 

J 

{ 

Vo ASSURE IRS 
J 

: 

{ 


$° 7, °CEL) TRARSTION MATRI2: 7 
$ ACS A4) 
FAR 


5] 


* 


Soma Gnneend bated 
* ; ‘ 


i 


HM eo 6 ke 


Subroutine READ 


HAHL «MO O4--B1 23 38:2 PARE OQ! 


SURCIT THE READE CHUM. HUMURP. ISTAPT.ISTOP. SRICT) 
COPAY URPBLE (3808) . EPROP «20S 

DIFERSION X29. 18) 

TTEGER URRELE SUB ICT 

ih i 

VEL CES '=4 

HM MPR 


MMe T=1,700 | 
CALL INPUT. HELOCY MRR. 2) 
IFCMPYEC.-2) Gil TO St 
(Gh Jette 
ICUS THEA | 
IFCISUB. UT ISTAPT) 10 18 
IFS }SLEA GT. IST0P? GA Tet 19 
HSIEH eg eps | 
NOPE E HIN SVE 1) 
ti COT HEE 
a ctemite 
“GUT TNE 
fee 
1 gfe =] A 
UPITEC?, 7098) 1 URPELESK) 
fe COTE 


MESS ROT TS TH 16) 


Het 


Goad 
* 


mea eed, 
ry ry ‘ a 


oe 


Subroutine INPUT 


FCP TRA WU HBIA-1 «HON @4-MAY-81 022382 15 PRE OF 
31 SUBROUTINE TNPUTCX, NBLOCK, MARK, ND 
er INTEGER. DATA. X 
ORS DIMENSION DATAT 258) ,X(29, 18) 
een4 1P=8 
ORES IPL=@ 
OOK 00 16 J=1,16 
ena? DO 1@ 1=1,29 
sd X(1,3)=8 
O89 = «16 CONTIN 
MIG HFRK=8 
Mall RERDCN, EXD=581 -ERR=681) DUI, TOUNZ, IDUNS, 1DURE, IDL, IDUNE, 
(DATA 1), 11,258) 
O12 HBLOCK=NBLOCK+1 
BIS K=@ 
€ MAIN LOOP 
mis DG 108 31,16 
mi K-K+2 
C STUDENT NUMBER 
@i6 X(1, =DATACK) 
C PERIO 
B17 KK42 
18 X(2, J=NODCDATACK) » 256) 
C SEGMENT HUMBER 
aS Oe) KK+18 
een 00 28 I=1,2 
21 K=K+} 
22 NIK=MODCDATACK) »2) 
23 X(3,J=X(3, J 4NTK 
924 DATACK=DATACK)-HIK 


25 =. 28 CONFINE 


K-12 


27 = LPSSDATAKK) “256+ 1P 


OZ IFCIPL.LT.LPS) GO TO 46 
MES IP=]P4+2% 
ARS | @ 10 B 
O32 46 IPL=LPS 
33 X(4, J=LP5 
C STICK X CROLL INPUT) SORLING [5 INCHES TIMES 168 MAK=6.24° 
OE4 00 3% 155.8 
635 K=K+i 
O36 Y=FLORT (DATACK) 196, 6198644 ISIGKS, DATAKK) 16.6 
0837 XL, J=INTCY) 
OCS 56 CONTIME 
C STICK ¥ (PITCH INPUT) SCALING IS INCHES TIMES 108 (O6. 335” 
39 DO 68 1=9,12 
GAG K=K+1 
0841 Y=FLOAT CORTAKK) 196. 6193234 ISIGH(S, DATACK) 16.6 
42 XL J=INTCY) 
0343s 68 CONTIME 
44 00 78 I=13,14 
45 KeKel 
C SIDE TASK START FLAG 
46 HIK=MODCDATAKK) »2) 
B47 DATACK)=DATACK)-HIK 


Rw pee Ee 


a ee 


ae] 
‘ 


Wee me 8. 


RO, 


ere 


nee 


Subroutine INPUT (Continued) 


* FORTE 1 HOIfri | PO O4-HRY-81 B22 58 15 PRE 682 


CHS 


* WR 


Gow ee 
. . 
oF a 
ot 
' ee 


X26, JD=EL26 J) OBNTK 
{ PEDAL POSTION SCALING 15 INCHES TIMES 168 MAK=3.25° 
VoFLOATCDATACE? 90, 0991+ 1SIGKS, DATACK) 18.6 
KO DDSINTy 
Th COTTE 
C PITCH INMCATED MAXSSG DEGREES. SCALED TIMES 188 
(OB [=$5. 15 
=4'4] 
CATH KOSDATACE AD CORTACK) » 2) 
NOEL CET (DATA EY 98, S491 SIGS. DATACK) 9716.8 
KL DYSHE 
a COTIME 
{POLL INDICATED [AMS DEGREES. SDALED TIMES 10% 
(G Mt y=IP 98 
vaiel 
TATA SAT ALE ALE DATACK) 2) 
VALOR DATE Hott, 498241 SIGNCS.DATACE 18.68 
TAD Seshteys 
“a COAT TLE 
1 HOTT TUDE 
Y=y'44 
eo DSDaTAY) 
Cee i Th feat PY 
Vaye} 
Wt. Si EY Se 


CHI PDSRED MAGEL4E. 5 FITS. SCALED TIMES 168 


a 
3 


Go, 
ined 


Pur 


fia) 


(re, 


srr 
ier 


4 


(4,"4 


Caed 
° ' 


(ay 
SAS 


BLZL DRE OCDATACKD 256) 


”. MEASING ft Ti 7A8 DEGREES. SCALED TIMES 108 


Vevey 
Hz. DY SFLORTCDATACE) 48, 549524 1S IGKCS, DRTACK) 918.8 


C VEPTICAL VELOCITY 


r=t'+} 

MZ DDATACK) 
{ 10 14.9 DEGREES. SOPLED TIMES 168 

ie Aa 

22.24 D=FLORT CDATACE 98, 045484 IS TGCS DATACK) 7108.8 
( CLARE DEVIATION 

Kael 

EXD3» J SFLOAT CDATACE) 198, S493 24 1SIGKS CATACH) 18.8 
{ S10E TASH MUMBEP 

K-21 

526. D826. J4NDOCDATACK . 256) 
( GEREPATC® RECYCLE TIME 

X23, D=ADATAH 26, 
( STATION SELECT TIME 

FKAZZ 

AZ?» D=DATAG/26 
C MIHE RESPONSE TIME 


(we X(2B, JAD CDATAKK) 256) 
_ OFS 10 COTINE 

[ wre HPK=16 
© Bie FETIP 


(es, 


ee 


Sil CATIME 


en PEUIID 2 | oe 


Subroutine INPUT (Concluded) 


, * EVE THIRD MON O-HAY-61 62238215 PRE Om 


wes, {12-2 
fara FETUS 
(5: O61 CONTINUE 
sr} PARK =2 
(S47 RETURN 
(i043, hb 


4 bed 


; 4 I 


38 


= 0 ees 


‘ rf Table 6 


PHAN3.FOR 


Return Cell Values for Flight Data 


ame? 
* 


| impose: The user enters error mean, error deviation, 
mean error rate, and the deviation of the error rate. The 
computer reads the data file and computes phase plane cells and 
outputs a count matrix. 


i 
Sones t 
* 


Computer Output User Input 


@ Run PHAN 3 


Same ¢ 


Ni Sail om ee Sas a as Se 
4 
Buen i 


<CR> 
; Pause -- Hit return and name <CR> 
. output file ‘ 
\ DK:TTEST1.DAT 
: ee Ps eee <CR> =a 
es Pause -- Insert data disc in drive 
a4 one and hit return CCR) 
j aie Pause -- Hit return and name <CRD 
ie Mie tnput file *DK:TTEST2.DAT 
oe Pe ON ae oA CRP oy J 
i i. Enter subject number 20935 
: i. <CR> 
‘i Enter check value to be analyzed: 
a airspeed — enter 90.00; altitude - 
2000.00 
j enter 2000.00; heading - enter 
- 90.00;rate of climb - enter 0.00 __ ____ SGR2__ ; 


‘ 
tae acct arene aM, DER rnc a te, cat ea 


. 
i Table 6 


j PHAN 3.F OR (Concluded) 
i : Return Cell Values for Flight Data 
{ Computer Output User Input 


TER ERMEAN, ERRDEV, 


RTMEAN, RTDEV! -340.56, 84.43, 0., 8.22 


ee ee ee een ee 
. 
i 
Baten i 
e 


R 
a or the variable to be analyzed: 
j pirspeed - enter 21; altitude - ‘9 
ie enter 19; heading - enter 22; 
a & ate of climb - enter 23 eek CR 
_ | 
, + 
<CR> | 


1, 9, -340.56 


CR) 


62, 73, 109.25 


<CR> ‘ 


INTER ERMEAN = Enter error mean 
RRDEV = Error Deviation 

RTMEAN = Mean Error Rate 

RTDEV__= Deviation of Error Rate = —. 


LP = Line Printer 


= Teletype 


Note: To end enter a negative number for the start position. 
You may enter more than one set of START, STOP, 
ERMEANS at atime. Remember to <CR> between 


40 


I 

1 

I 

i 
rob, 
t 

t 

c 

“e 


Example Output PHAN 3 


QRIFCT: 28S 

VR IARLES 19 

WEEK (F SAOPLES: 9 

VARJAELE SCALE: «=, ABA 

START PORT TION: J 

SIU PORTION: 9 

PEA (¢ FRROR = -M40. 5688 

HER DEVIATION = 4. QR 

IEA (t ERKC® RATE = 0.6088 

DFU ATION OF ERROR RATE = 8.2208 


CHI URLOES 


8 wee - 
2 


ERE HE 


\i 
iI 
ie 
I 
I 


eer eee eos 


cn] 
2 


ft. 

URIE(TS 7A 

WRIHBES $9 

HARTY (¢ SLES: 

UHOHAE SCHLES }. 00888 

Sie ASITION = 62 

Sk ASIC = (78 

HS (¢ FAT = 109. 2388 
FRR DEAT ION = 04.48 
YEARS ( ERROR RATE = 0.0088 
“TELDAT ICN CF ERROR RATE = 8,708 


7- 


pag ponte pieceere 
e . . 


Pe) 
als 


Shi) SALES 


Example Output PHAN 3 (Continued) 


cee va ee RE ee oe == ee 


a a rr 
: ? 3 4 5 


Ah -8 «45° 41.5 


J -39,258 12,008 6 
7 2.750 14.000} 
% -13.750 16,0@a 7 
4 2.78 90? 
oNS,78@ 0A 1? 
6 5,758 4,008 12 
7 7,780 -4,008 12 
& 378 ~7,00R 
4 1.750 3.0081 

1M 4,758, 40m 7 

1) $3,758 8,07 


UIE: AAS 

UMP IAF: 99 

MMER CE SAPPLES: = 28 
UHRIOPLF SCALES = Ree 
SUR) PUSDT TCR: rt 
oe MTTINS = JAR 


HEA (CF EPRTR = HBR 

FRET DFAT ION = 94, 4308 

NEA ( FRO PATF = 0, 0RER 
DEAT ION OF FRROR RATE = 8.2288 


CH URCUFS 


eee OPT DESY Carey EEE 
j ? 3 4 3 
Cee Ae eee 
6 7 g 9 «8 
Se ee en eee 
" 1? f3 4 «5 
Se ee 
tr, 17 §8 19 ma 
ST a 
ae a me 74 
ee ae er See a 
Pac -5 #5 448 
FRAY 


w, 


42 


a tc a es 


ok A a a a a mc 


(a 
& 


Example Output PHAN 3 (Continued) 


e [ SHES ERROR RATE CE ERROR CELL 


| 1 $8,118 ee 19 
a I, ” $37,518 R19 
3 $79,110 4.08 «14 
c 4 SA16 2008 19 
| i 5 N68 -10.0@ 19 
| eB 16118 em 19 
| | ytd Me yey 
- Ne “ae 19 
| i 4 99,110 -71,0 
: "se 67.110 14,008 DB 
. 13 43,418 rym 
J. se ee 
a> 7 1% -7, 28 71.008 R 
14° -77.890 1.0 oD 
th 44990-1600 
i aes “11,008 1? 
| Ve 71.8% “12.0817 
' 1@ 60,898 “15,008 2 
14 9, 8m “14,008 71 
S- " -12,808 RO 6 
: 1-1, 8 em OH 
: oP S57, 898 2,908 é : 
| 8. ASG, 11.900 é , 
; “HOR ha, : 
i * 5, 8 1%, 0G 2 
f BAe BOR nme 2 
, ae Ss: ime 


AMIEL AOS! 
MESHES J 
MMEEP (t SAMPLES: 7 


eT o. 


; METHME SCHLES 100888 | 
f Sib) MEITIOMS 192 
H { SOP EUSTON 12% 
MHS ( FREOR = 96.4288 
3 Het HIMATION = 2, 408 
: I MCh 1t EMA® RATE = e.eane 
TEATS CE ERROR RATE = 8.208 
E CRY IES 


i 


Ganeg 
. e 
> 
ee) 
> 
as 


Zi = Se ae eet ee ee RR 41.5 
Paes . ix 7 g 9 18 


eur 
i} ¢ 
' 
— 
4 oD 
t 
t. 
= 
~ 
ao 
-- 
wo 
t 
_ 
Py 
mn 


“3S 2 BMD 43 


ar 
| 


-27,.428 

7 77.48 
o 75.4 
“$4.48 

-7 4 
-4,4.8 
#.88 

Am | 
16.390 
19,0 
“158 


fl 
[ 
i. 
i. 
i. 
{. 
I. 
I. 
ie 
I 


I CHL URC LES 


7.08 
ee 
4.0m 
4.088 
5.688 
7.0 
9.08 
3.088 
1.08 
1.906 


% ERMOR RATE CF ERROR CHL 
J 


~ 


ee 
NNN DNR 


Sos 


Example Output PHAN 3 (Continued) 


Example Output PHAN 3 (Concluded) 


SOOTHES GOSEGEGTEoeGooGoqnoaenns 
OeeEeeaaeaeceoeeeaetoaoreuaae s 
Sear Gee aecaueceqaeaearocta ean arn & 
SOOQGeT aGoee qe ea eeaGnatoqanr- & @& 
SOSGSeeEe oe ureaeaaeeeaeoeroeraa-ae ae 
Eee ceqecreeeqeeeeaanraoucaqneaana 
SSQEer aa eeeoueanaqaecaanra aaa cae 
SGSOEGe te eee anaaceaetee aac 
SSG ooo ereoaqaaeaenaaeqeaqqarqe ae 
SSSSHSESSETESEELEOEOSETEH SOEoE 
S@SGESGHOOEe Ge eGawreaceqqoaceq 
S@SGSOSTO@EqEGQoreeaearaaceaq 
SOG etree GOeooeenacaeacteoa 
SSESSSESMESSESKOSEGESETOGaGESeES 
S@Eeoeoeaeoeeeeeooeoeoeen-Getaet eae 
Se@Seceoeeteeq@ OOo aqqecoeceacaeace 
See eoe@eeteeeooeoee2rotaaccecea 
S@EeeeeceeGGeoeoeoneceereegcqceraca 
fe-GQeeevrvee@eneoeoeeatecaoqecece 
@SSeceo- SF ESTOE--Seseeeereoeact onan 
SBOE ER ESSTSSOTEESSSESE SE OSes 
Pd meee eae wen ee seus 
Se ee ee SESE SES ESS SOGCEESEESESES 
NNO SSESEHEEGCHETEBSHESESCEWVSERQeEaESCAESE =| 


LSESSNSESETESSESESIFSEDTSCSGS SSS 


SLESTLERELEEEEREeeeses 


SSSoesesesaoade SHSesecsesd as sss 


SCSETESTSESSESLSSSS SELEE 


s es s s s s a 
SOS aoe ooaqcorFequceaqoaqcoacoua 


SSSSSSRSLE SESE eSRResErce 
SOEeCeqodoesaocoOooquwucoaoeqocac 
SRLELLSLE HS SSE RLS REREKES 


s s iJ = s a * 
SOeoQVeOe Gago Foaaocaoreaoucaae 


Soo gan ruaseeseoeqe ase coooaan 
S@EERESEKESRELEREEESE SESE 
SESescseasasssodcisgcEescsiscedid 
SEEEREEEREREECESRESERESEES 
GESeaedrasIasaeseosc si iaicGcasssisd 
SSSSESEBEKREEERERERESKBEE 
SSEeeesseoeseiaceaqicgasssicscsaavad 
SEKRSEEEKEREREEEEEEREEES 
SEH Sdsrscasceciseicace awe 
ESEESEEEEREREEEEEEEER EE EE 
See eaceraasaecacaccide ado d ad 
SESSESEEHHREEESEREEEMEREES 
Sessa aedidsec cdi deesada cel 
SHCEKEK ELH ERER ES HS EEE EE 
Seceseeeeeaeialaeaearusee he cc 
EF ELUTVERTCRLE HH SSSEEE HE 
SGSGSGSEEeEeEiasecsacaasecaceaacd 
SBFHKESSESEGEEEKEKEEKEEGEEES 
Ges ecee seceded didccecde adc 
PESSSKSSSSESSHEREKEEREEEEE 
ESGHesecseecce dence dackc&ceEaccede 
TECSTFECKCKLEST CHTEESSTLES 
SEs eee Gece s asec acedc 
SSSESSSEESCESSESEREESEECEEE 
SSEser sre esse sites sid se SSecocea 
IRESSEFESSSERSESETESEESESE 
ESESSESSSSSESO SSS SSeS SESS 
 FEKEMESSTESEESCCSKESELES 
Seer seece een Se qaqaeaeqaeaaaa «ac 
SESC ESS TCCSSSRE KE CEKE SE EE 
Sees aeeqoaeeeGe scence ecaacecaec 
SEFSTESESTSESSESTSEESTEESSEEE 
SSE t Sees ses acseeecaccece 
SET SSE SSCS SEES SSE SERSEEREE 
SSSSESSSSESSST SC SSeS SSeS acces 
ESSEC SESE SHSETESET ESET ESELRE 
: Serscececeecrtsse eeeeacqgace 
PESEGESSSESEESEEEEEEEEEEE 
Sosrcrsssrewvs 2 Perce d 


<enngpe <- aps 


Buociatn. § Gens t. 
, i 


Gu 


ELELLELL 


Beene 
- 
-_- 


‘ 
ce 


Program Listing PHAN 4 
FORTRAN 1V Wife} | WE ZL Si GG 1528 


PHANG. FOR 


PURPOSE TO READ MULTIFLE COUT ATRILES E- 

TO FORM CE TRAE TTION MATRIS 

DIMENSION KIMAIEC. AD DIMRE 8 AMT) 
065 1=1.0 

06 JEL 


> i 
o 
oo 
a 
m 


PARE HiT RETURN € HOE CTH FILES 
CALL AGS IGN. SoD NER NO hd 

$1 PRUGEH HET GPE THA FE 
CALL AGSIGEZ. ERNE ds 
0G 16 1=1,4 


00 13 J=1.4 

KTATXC DSA De: 
13 OONTINE 
12 CONTINE 
C COMPUTE THANG TICK NATED. 

(G $4 121.4 

SUK. 

MH IeLA 

ATMA D=ELCATPAT 
1S SUREQMHATIATS Tt 

O15 FLA 

JECREC.A) SURE, 
Ye AIMATACL DeATNAEG 


14 CONTIME 
PAUSE HIT PET SHELL SFG COTM RGR Ste’ 
REPOCS, Lol TAS 


CALL CUE tZ 
TFCIARG. ECT G1 
it Van bd eee 
$6 WRATEC LAIMA Ds 232) 
CALL CL0GECI) 
STOP 
$@ FORT C14, 
10) FORMAT CII: 
1@ FORA CA6S. 2: 
END 


| 


Be 


Go nsce § 


¥ 


oe Be oer oe El oer lll one 
e 8 @ @ «© «# 8 ¢ 


Table 7 


PHAN 4.FOR 


Purpose: Reads multiple count matrices 
and forms one transition matrix 


Computer Output User Input 


@® Run PHAN 4 


<CR> 

Pause - Hit RETURN and <CR> 

name Output File *DK:H491T3. (IAT 
CR> 

Pause - Hit RETURN and <CR> 

name Input File *IK:H491C3.DAT 
<CR? 

Pause - Hit RETURN and <CR> 

Input 1 for Continue, O for * 1 


Pause - Hit RETURN and 
name Input File 


«CRY 
*DK:H201C3. DAT 
<CR> 


Pause - Hit RETURN and 
Input 1 for Continue, O for 


Ne a nee ee em en eons, 
cy 


- FURTRAN TU 


PO 


7 


Program Listing PHAN 5 


HBIA-1 | THE 28-UL-81 1:06:54 
PHRE. FOR 


PURPOSE TQ READ MULTIPLE COUNT MATRICES ABD 
TO FORM OME TRANSITION MATRIX 
DIPENSICN KTMATX(ZS, 25>, TTHATEKZ. A) 
005 I=.A 

M6 ELA 

KTMATXC), D=@ 

CONTINLE 

CONTINUE 

PRUEE’HIT RETURN & NAME CUTPLIT FILE” 
CALL FESIGHCL,“ 7yc1, NEW NOL) 
PRUSE“HIT RETURN @ NAME INPUT FILE’ 
CALL ASSIGK2,° *.-1, ROO". HE’) 
D0 18 }=1,2 

REPOCZ, 1) CIMA) Ish) 
0012 I=L. 

06 13 Jz, 

KTPATXC 1, SRTMATEC J+] HATA DD 
CONTINLE 

CONTINLE 

PRUSE"HIT REVUE SIKPUT FOR CONTIIRE.@ FOR STOP’ 
REAOCS, S00 TS 

CAL C.05EC2) 

JFCIANS.EG.1) Gh 1G Mt 

DO 17 3-124 

WRITECS SOCK TATU, J=1,29) 

CALL CLOKECI® 

STOP 


18 FORMAT C14) 


4 18 


FOBT C1) 
EN 


PRE 081 


Table 8 


PHAN 5.FOR 


Purpose: Takes multiple count files and creates a 
single count file 


Computer Output User Input 


| @ Run PHAN 5 


€R) 


Pause - Hit RETURN and <CR> 
name Output File *OK:HPC5. DAT : 
<CR> 
wees - --— + — we eee ee | 
Pause — Hit -RETURN and <CR> 
name Input File *DK:H991C3. DAT 


<CR> 


Pause —- Hit RETURN and 


<CR) 
Input 1 for Continue, O 1 
for Stop <CR> 


Pause - Hit RETURN and 
name Input File 


CRD 
*DK:H201C38.DAT 
CR 


Pause - Hit RETURN and 
Input 1 for Continue, O for * re) 
Stop <R> 


Program Listing WTMAT 


~ om 


> IRAN IU | HBIARL «TRE 20-RR-CL are9:51 PAE et 


t ( PROGRAM ‘WIMAT. FOR’ 
a C CREATED FOR THE HLIA SYSTEM 
( ORIGINAL DATE: 4~JAN-90 
= C PURPOSE: STOCHASTIC ADJUSTMENT FOR THE ‘Me’ SISTER 
La C TRANSITION MATRICES AND THE SCORE ASSOCIATE WITH 
C EACH MATRIY: ARE THE INPUT URPIABLES. 


= C OUTPUT IS A WEIGHT MATRIY WHICH MAY BE LISED Tu LCE 
; C THE PERFURMPNCE (F A SARPLE. 
= (al CORON “PEC” TRARS. TSPROR, DLIMIT. 
i DLHCL: WEIGHT, 
? Z PACORE .COLCHE . GAIN, 
a 3 FC TINT EPRACI PPL IM. 
4 ELOG!.WHGH 
we? COMO TELQCKS FRITS. FITZ. ERTS. 
: i ELF MUMAT I TMGIE, 
2 MPODL] -MARALT  IALGHE.. 
3 MAINT 03h). LENGTH 
4 TRIG! JHIGE. FEY. 17 
= OR OTPENSIOK TRAE CIa8) : TE S35: : 
1 DUMITCGZS) LHL e cS) 
. 2 IG EA:, 
wb 3 PSCUREC TS), CSCORE Ch SAI MLPA), 
4 IFAT ERTS . 


( READ CONTROL VARIABLES 


I te PRUSE “INSERT DATA G13. FULMER ic st HIT REVERIE 
i CALL CURERD 
C INTIALIZE WORK UARTAELS 7 
(He. 


] he CALL INTIAL 
(READ THE PROGLEN URRIAELES 
wer CALL PUREAD 
( COMPUTE THE LIMITING DISTRIETIUNS ECE 
] ( THE TRANSITION MATRICES 
eee: CALL THELIM 


C COMPUTE THE THAISTATE PROERELITY HGziLE! 
I 9ee9 CALL TPM 
CCUMPUTE THE GAIN FACTOR FOR EACH HASRI 


mia CALL TMGAIH 
C PERFORM THE STOCHASTIC ACJUSTNENT PROLELUPE 
ae @1! CALL ADJUST 
C OUTPUT THE FESULI: 
~ 2 CALL ARITE 
GIS STOP “END CF PROGHHT STCHDI. F(R’ 


EO 


I . Subroutine CVREAD 
_ FOO THEE TU RIAL = TUE GS-RPN-81 03258257 PRE Ot 
eo at SUEPOLTINE CURERD 
( PURPOSE 15 10 READ THE CONTROL URRIABLES 
> een COMMON ARELOCK/ TRANS, TSPROB. DLIMIT. 
j : i DUHOLD WEIGHT, 
2 PSCORE, CSCORE. GAIN, 
Z 3 POJINT ERPADS -ERRLIM, 
| 4 WL.OU, WHIGH 
* OARS COON “TRLOLY IFRETI. IPRHT2: IFRET2. 
! TOLFEDWPMAT . ITHSIZ, 
; 2 HRODLLHRORIT, IRDOE, 
i 3 HSIET, ISOM, LENGTH, 
4 I81G! . IBIG2. IBIGS, IPSMT! 
(ness DIFENSICE TRANS(2202) . TSPROB( 2088) . | 
1 DLINIT408), DLHOLD<488) , 
2 VEIGHT (490), 
3 PSDTRE (28), CSCORE( 28), GAINC 28) . IOLFHD( 28), i 
4 IFPITT 128), IFRAT2(28), IFRMTS<28) 
- + WS WITEC?. 7189) 
Man, (FEL ESIGN. 6-1, °R00. WEDD 
(97 WPITE’?. 7118) | 
won (0H CL ABSIGK2.° “o-1. HEV “NED | 
(ses LTTE? 712) 
t (sts: PEACCS.S1GQ) HMMAT 
oe il WITEC?.7138) 
(at PEAS. 5108) ITNSIZ 
_. MS WITEC7.7142) 
tats PERO S 5209) ADSINT 
ee ts, URITEC7. 7150) 
W31e. READS, 5288) ERRLIM ; 
7" ay WITE(7, 7168) 
we HIE PEADC5, 5208) ERRADS 
sits WITEC7,7178) 
mo C2 REPOS, 5108) MRODLI 
S wazs WRITE(7,7198) 
** wiz READC5. 5108) MAXAIT 
0B WITE(7,7198) 
i 4 READ(5,5308) CIFRMTICE), I=1,28) 
. 6B WRITEC7,7208) 
126 READS 5308) CIPRMT2C1), 131,28) 
wu? WITE(7,7210) 
| woe READS 5308) CIPRATSC1), 121,28) 
wz WITEC?, 7228) 
mes PERO 5288) WOW 
| OGL. ss WITEC7,7238) 
ez FERO(5.5208) WHIGH 
ec; WITE<7,7246) 
| er PEAS 51G@) IPR! 
oo PETLEA! ‘ 
5108! FORMAT C17) 
6k: FOPURTCFI3.6) 51 
60% 725 FORT’ WARE THE INPUT FILE’) 
SESS FORME’ HRP THE RESULT FILE’) 
ie R41 Th FOMATC’ ENTER THE NUMBER OF TRANSITION MATRICES (1)7) 


[ Subroutine CVREAD (Concluded) 


ROC ER YL! HeIAR-L = TE GS-81 63:58:57 PRE 922 


J ese 7138 FORMAT” EATER SIZE OF ROW OR ONLURH FOR MATRICES <1)“ 

C3 7140 FORARTC” ENTER THTIAL UPLUE OF ENTRIES IH GEIGHT MRTRIX (R)”) 
_ 44 715K FORPRTTC” ENTER ALLOMFBLE ERROR IM LIMIT COPPUTATION (R)’) 
[ @85: 7toe Feet” ere pLLOUELE ERROR IN SCOPE DEVIATION @R)”» 

- (686 7170 FORMAT” ENTER PERIPLE INTERATIONS FOR LIMIT COMPUTATION 1D’) 

COM7 7188 FORMAT’ ENTER PRIMER INTERATIONS FOP WEIGHT AMRUSTHENT (1)°) 
| GENS 719 FORPTC” ENTER FORMRT FOR TPAISITIOM MATRICES READ (20R2)°) 
{cosa 708 FORURTC’ EATER FORYET FOR SCORE READ) <20R2)”) 

cent: THN FOPRAT(’ ENTER FORPRT FOP LEIGHT OTPUT (2892)°) 

Gato 728 FORMAT’ ENTER LOWEP LIPIT FOR GEIGHT PATRIX CALL (RD 
720 FORMAT” ENTER HIGH LINIT FOR WEIGHT PATPIX CELL (RY) 
"GES 724g FORTATC’ ENTEP INTEPATION PRINT SUITCH FOR PORISTMENT (1)’) 
gee BD 


boa 
=e 
nh 
ts 


1 
| 
| 


4 hava 


6s 


oe ee ee 


Sineg eaeeg Said 
e % ® , 2 ' 


se 


Subroutine INTIAL 


F 
we 

_ RORTHAR TY HBIA-{ = THE @S-+HFAV-€1 08:51:51 PRE 68! 
i ABEL SUEROUTINE INTIAL 


( PURPOSE 1S 10 INTIALIZE URRIABLES IH COMMON 
62 COMMON /ARBLOCK” TRANS, TSPROB, DLIMIT, 


! 1 OLHOLD. WEIGHT, 
She iiss #28 PSCORE-CSODRE GAIN, 
3 POSINT ERRADS, ERRLIM, 
: 4 WU. HIGH 
+. 3 COMMON /TBLOCK’ IFRMT1. 1FRAT2. IFRATS. 
1 1DLFHD-HAPMAT, ITHEIZ: 
2 Men -PRERIT, IADGHE. : 
7 3 HAINT, ISOTH-LENGTH, 
4 TRIG - IBIG2, IBIGS. IPSWT1 
A DIMEMSICE TRAGS( 2008 . TSPROB( 2088), 
1 DUINTT¢402) .OLHOLD(400) « 
2 LE1GHT(490), 
z FETE S22). CSCTREC 20) .GATN(28), TOLFNDC28), 
4 TFPOT 128) TRAIT 2628). IFRMT S22) 
-- ON IBIG!=20"0 
(eee, IBIGZ=408 
+ UA? IBIGI=28 
woe 8 IPOHER8 
weirs HAINT=@ 
eats ISOTMELTMSIZ#1TMSIZ 
ttt t LENGTH=NUMMAT#1 TMS 1281 1512 
Hz OLINT=1. @/FLOATCITMSIZ) 
at” D0 108 I=1, 1BIG! 
mitd TPRHS(1)=8.8 
eats TPES 1)=8.8 
sts, IFCL.GT.IBIGZ) GI TO 198 
ss} VEIGHT CI =ADJINT 
aK DLHOLDCI=DLINT 
had DLIMIT(1)=8.6 
9 | IFC1.GT.1BIG3) GD TO 108 
a PSCORE(1)=@.8 
24 CSOORE(1)=8.8 
5 Fas CAINC1)=8.6 
we. IOLFNDC1)=@ 
27 ~=—«108 CONTIMEE 
erze RETURN 
irs 1) 


I 
I 
c 
2 
; 


[ , Subroutine PVREAD 
. 


ECETHA IU HRIREL «TUE GS-81 @3851227 PRE 801 


j oft SUBROUTINE PUREAD 
C PURPOSE IS 10 READ THE TRANSITION MATRICES AND THE 


: ( THE SOORE AGSOCIATED WITH EACH MATRIX 
| (GET © CORPICN /RELOCK/ TRANS, TSPROB DLIMIT, 
: i DUNOLD WEIGHT, 
2 PSCORE, CSCORE, GAIN, 
{ 3 ADSI ,ERRADS ERPLIN, 
e 4 WLO!, WHIGH 
ES COMMON, ATELOCK IFPRTL, IPRHT2, IFRITS, 
, 1 IDLFHD, MUPMAT, ITHSIZ; 
| ae Zz HAXDCL_LNRERIT, IADGHE, 
3 HAINT, ISQTH, LENGTH, 
4 IBIG!, IBIGZ, IBIGS, IPSWTt 
i 04 DIMENSICH TRAES( 2008), TSPROBC 2808) . 
1 DLINIT<408), DLHOLD(468), 
2 VEIGHT(408), 
| 3 PSCORE( 20) CSORE (28), GAIH(28), INLFHDC28), 
: 4 TFREIT 1 (28). TFRIT2¢28), TFRMTS(28) 
085 IK=MEMRT#ITISIZ 
2° 8S Ist 
owe 00 28 KEL IK 
eae VET +1TH612-1 
er PERC  IPRET 1 CTPARSCL) »T=12, 12) 
eniics H=1i4ITHS1Z 
i! 2 COMTIME 
tz PEA, TFRAT2) (PSOORECI, Jef HUPRAT) 
eats RETURN 
ais Ew 


& pas : Pe See a 5 S% Re 
rege eae rca i aa are 
Ano ee ere = en aie ae A ta af Pa  T E 


Subroutine TMDLIM 


|. py 


. FOOTED TY HRIPRE TE OSA ORS PAE OM 
j » eee LERCLITINE TROLIM 
C PURPOSE IS TO COMPUTE THE LIMITING DISTRIBUTION 
; ¢ FOR EACH CF THE TRANSITION MATRICES 
j tz COPMON “RELOCK’ TRANS. TSPROB,OLIMIT. 
1 DLMOLD, KEIGHT, 
3 2 PSDIRE  CSDIRE, GAIN, 
| x ADTINT ERPADY ERPLIM. 
-- 4 ULC! AHIGH 
OS COPMON /IELOLY’ IFRATI. IFRNT2, IFRMTS, 
| i TOLFID.MUPMAT, ITHSIZ; 
st 2 MRODLL MRRRIT, IRDONE, 
3 NRINT, ISTH, LENGTH, 
4 IBIGI, 1B1G2, IBIGS, IPOMTI 
| ee DIMENSION TRANS( 2083) . TSPRNRC 2008). 
1 DUINTTC408) ,OLHOLDC 488) , 
z 2 EIGHT 498), 
{ 3 PSOLPEY 20). CSCORE (20), GATN(28), INLFND( 28), 
a 4 TFPNT1¢ 23), IPRMT2(28), IFPMTS(28) . 
§ THEO 
1 SHS NC 568 I=f OTL 
Qe ts; my Tee tf AZER r 
Mikes IFIDLAY D308) G10 me 
wots VSTAPT=CT~{ HT TMSIZ4 
(ance PSTOPSYSTRET+1 TMSIZ-1 
est CEPR Hee 
wat? 0 208 R=AYSTAPT.RSTOP 
He LETAPT=CK=1 41TH IZ# 
(std LETCPSLSTAPT+ITIGIZ~1 
feats 00 103 L=LSTART.LSTOP 
(ats THOLDCL)=0LHOLDXK)#TRANG(L) 
wit? WRITEC7.12) THHOLOCL),OLHOLDCK), TRRRE(L) 


M2 12 POPMATCSFIG.3) 

59 $08 CONTITE 

fs 20 OUT IE 

324 KOR 

Lary 00 228 K=KSTRPT STOP 

O23 KDUREKDUH I 

es24 LSTART=(J-1 9 1SOTHHDUM 

a5 LSTOP=_START+1SOTH-1 

CHS DLIMIT(K)=8.8 

WZ7 00 216 LaLSTART.LSTOP, ITH5IZ 

OL DLIMIT CKD. INITCK)4THHOLD(L) 

79 WRITE(7, 13)DLIMIT(K) 

36 13 FORPAT(F 12.4) 

O31 = 216 CONTINE 

O3z ERROR@OL IMT T<K)-DLHOLDCK) 
SSRERR=SSUERR ERROR VERROR 


“| 1 ot 


DLHOLD(K)=DL IMITCK) 55 
C5 20 ONNTINE 

S% IF(SSOEPR.LE.ERRLIM) 1OLFHOCS)=1 
22 380 ONTIME 
ong IFS 
aes 00 400 MeL HUPAT 
IFCIDUFHDCH).E.@) IFSUMELFSUM 


om f.. 
32 


bo 


B. 


. FOR TRRE 


HOIR-1 «= TUE GS-01 51:48 PRE 972 
408 CONTIN 
IFCIFSULEO.MRPRT> GD TO 988 
388 CONTIN 
WRITEC7, 7108) 
O 10 999 
403 CONTINUE 
WRITEC7, 7116) C1, TOUPNDCT), I=] -NPPRT) 
WRITEC?7,7128)  CS.OLIMITCS), J=1. 150TH) 
99 CONTIME 
RETURN 
713 FORMATC’ FAILURE TO FIND ALL LIMITING DISTRIBUTIONS’ 
i] ” PROGRAM EXECUTION CONTINUES REGARDLESS’) 
7126 FORMAT (’ SUCCESSFUL COMPUTATION OF LIMITING DISTRIBIT IONS’ ¢ 
| * MATRIX’ 3X.“ TERATIONS” 
2 QAC/IK: LTR, SY) 
7128 FORMAT’ VECTOR FO® LIMITS FOLLOMS: “- 
1 TOES 6. LIMIT’ 
1 We ES TS. 2K,F9,9)) 
BD 


Subroutine TMDLIM (Concluded) 


7 Re EO - 


ent 


Qe 


ws 
"(Eo 
i we 
i G63 
—p (eS 
> eetit 
at 
en 
~ ett 
aw Mid 
res 
"we fete 
4, a7 
rene 
ae 9 
_ $e 
t 
[ e272 
| ears 
Ps 
a dd 


Subroutine TSPM 


“ROR IU HIRE TRE GS-81 60251257 PRE Ot 


SUBROUTINE TSPH 
C PURPOSE 3S 10 COMPUTE THE TRANSTATE PROBABLITY MATRICES 
COMMON ARELOCK’ TRANS, TSPROB. DLIMIT, 
i OLHOLD, WEIGHT, 
2 PSCOPE -CSODRE, GAIN. 
3 AO SINT ERRAD. ERPLIM, 
4 SLO WHIGH 
COMIN “IBLOCK IFRATL. IFRATZ, IFRATS. 
| TOLFID: HEPAT. IPG IZ. 
v4 MAD -MRORIT. SADOME - 
3 HATH. ISOTM. LENGTH, 
4 IBIG! - IBIG2: IBIGS, IPT! 


DIMENSIONS TRAS< 2522) . TPPOB( 208A) . 
OUIMIT S400). HOLD 480) . 
HEIGHT 462), 
PETES 28) -COONPE( 28) .GAINC28) . TNUFHDC 28). 
TFPRT C28). TERMT2( 28). IFPMT C28) . 
STPPPEK 28) 
0 SG3 [=f EPNAT 
STPPPECT=8,8 - 
ISTAPT="1-1 HT TMSIZ+1 
SOTQP=ISTART+I NGI Z-1 
10 260 JEISTART, ISTOP 
ESTRPT=CI-f eT THI Z+4 
ISTOP=ESTART+1 HS 1Z-1 
0G 168 K=KSTART.KSTOP 
TSPROBCK=DLIMIT(J)#TRANS(K) 
STPPRECT)=STRPRBC 1)+TSPROBCK) 
123 CONTIME 
200 CONT RE 
260 CONTINUE 
WRITEC7, 7108) (STRPPBCI), 1=1, HAPRT) 
WRITEC7,7116) (1. TSPROBCT), T=1 | LENGTH) 


CS Bed bo 


RETURN 
7108 FORMAT(’ CHECK SUM FOR TRANSTATE MATRICES’, 
! ” SHOULD BE EQUAL TO 1.@ (APPROXIMATELY) ’ 
2 1X, ¥18.5/)) 
7118 FORMAT(’ TRANSTATE VECTOR FOLLOWS: “7 
i * INDEX’, 9K, “URLUE’, 
2 BOC/AAK, TS 1X,F13.5)) 
EN) 


Subroutine TMGAIN 


_ FORT IU = Hetret «TRE @S-HPV-B1 325271! PRE ent 
S. oe! SUBROUTINE. THOAIK 
C PURPOSE 1S TO COMPUTE THE GAIN FACTOR FOR EACH 
CCF THE TRANSITION MATRICES 
1) eae COMO /PELOOK/ TRANS, TSPROB.DLIMIT, 
1 DLHILD. WEIGHT, 
7 2 PSODPE .CSORRE GAIN, 
| 3 POJINT EPRADS  EPRLIM. 
ares 4 ALON. HIGH 
wer CORTON “IBLOCK.” IFPMT I IFRMT2. IFPHTS. 
1 TOLFHD HUOPAT, ITS. 
: 2 FECDUL PRRRIT. IADORE:. 
3 HINT, ISOTH. LENGTH. 
4 IBIG!  IBIG2. IBIGS. IPSHT! 
{eons DIMENSION TRAS( 2802), TSPPORC 208A), 
1 DIMI T<402) LHL 49) 
2 UEIGHTC400), 
3 PSNNPE< 28). CS00REC28) .GAINC28) . TOLFHEY 28). 
: 4 TFPRT 1 (28). IPPHT 2028). TFRMT3C2R) 
oes NC 200 T=1.4LPeAT 
BES, SCTAPT=(1-1 #1SOTHE 
ow OF JETO@= ISTRPT+ISOTH-1 
(S22 (160 JSISTAPT. ISTOP 
eee CAT =GAINCD)47SPPOB(.D#TSPROBC.DD 


mates 10 COTINUE 
Pat] 20S COTE 
Ld WRITEC7,7108) (GAINCT). J=1-HUMMAT) 
ohh PETUPK 
C214 7108 FORMAT(’ COMPUTED GAIN FOLLOUS:’/ 
1 41%, ¥F18.5)) 
ay 


Subroutine ADJUST 


HBIARL «= TUE @HARY-61 05:52:25 PRE Ol 


~ Si} SUBROLTINE ADJUST 
( PUSPOGE 1S 10 PERFORM THE WEIGHT MATRIX ADJUSTPENT 
{ PROCEDURE UNTIL CONUERGAICE OR MAXIMUM 
( INTERATIONS ARE PERFORMED 
re COON APELOLY TRANS, TSPROB. DLIMIT. 
OLHILD- WEIGHT, 
PSCORE CSOIRE GAIN, 
FOJINT ERPAD, ERRLIM, 
WLC! HIGH 
TFPNT 1. 1FRAT2Z. IFRATS. 
TOU. WOMAT, ITMSIZ, 
HFEDLI .MRERIT. IRDOHE. 
HAIHT. ISOTM LENGTH, 
TRIGI - IBIG2. IBIG3. IPST! 
DIMENSION TRA 2808). TSPROBC 2888) , 
CLINT (408). OLHILDC 4AM). 
EIGHT 408°. 
PEELPES 88). SSOTRE C28) | GRINC 28) » TOLFHDC28) » 
4 TERT TR. SEPA 2C28). TERAT C28) 
COE al (ORNS 
A bs co 
SAG T=] MERIT 
NATHTSHALHT +4 
00 38 J=1 NPAT 
CSODRE(J)=8.6 
DO 108 K=1. 1S0™ 
L=(J-{ 1 SQM 
CSOORE (J =CSOORE (JE IGHT (K)#TSPRORCL) 
106 COUT IME 
ERROR=PSCORE (J)-CSOIRECT) 
AUGAINFEPROR/CAING) 
00 208 MEL. ISATH 
HEC J~Le1SOTM 
WEIGHT (M)=RE IGHTCM) +AGAIHOTSPROBCK) 
TFCETGHT(A) .LT.ALOM) VEIGHTCA)=MLO 
IF CWETGHT <M). GT. UHIGH) VETGHT(M)=HIGH 
2G CONTIMLE 
208 COMTINLE 
SSOERR=3.6 
00 328 J=1 .NPRAT 
CSOURE(I)=B.8 
00 316 K=1.1S6T 
=(J-1)s1SQTHK 
CSODRE (J)=CSO0RE( J) 4NE IGHT(K)#TSPROBCL) 
U@ QUITINE 
ERROR=PSOORE (J) -CSDIRE(J) 
SSREPR=SSUERRHEPRORVERRIR, 


ee en > emcee a ome 
a 


328 CONTINUE 
IF(SSHERR.LE.EPRADJ) IADONE=1 
IF CIADOHE.E0.1) @ 10 3 
CINSSE=FBS<SSEHLD-SSPEPR) 
IFCCINSSE.LT.CLIMIT) GO TO 38 
SSEHLD=SS0EPP 
JFCHEOCHRINT . [PINT 1). HE.@) GD 10 Sa 


er 


4 


Subroutine ADJUST (Concluded) 


FORTRAN I! HOLA] «TUE @S-RY-81 43:52:25 PAE @02 
J -coae 338 CONTINLE 
$7 WITE(7,7108) NAINT, SSOERR 
~ ae DO 408 J1=1, HURT 
4 e649 URITECT, 715) J1,PSOORECIID, CSDDREGIID 
(re: «© 408 CONTIME 
- EL —-WRITEC7.,711@) 
i W8z 00 410 Ki=1, 1ST 
«+ FS WITEC7,7115) Ki VEIGHTCKI) 
(54 418 CONTINUE 
j cat IF CIADOME.E@. 1) GD TO cae 
£ . aei7 IFCCINSSE.LT.CLIMIT) GD TO 788 
65% SAB ONTINE 
EA WRITE(7,7128) 
i} BAL @ TO om 
"(062 688 CONTI 
n3 WITEC7. 713% 
(i:4 @ 10 988 
+ - O65 708 CONTIN 
A, WRITE(7,7148) 
"(67 OBR CONTIME 
<. me RETURN 
@6% = 7100 FORMAT(’ ADRIST PRSS: “.14,° SUM OF SQUARED DEVIATION = “.F13.5/ 
oe 1 ” INDEX’ ,8X, “SAMPLE? .6X, ‘COMPUTED’ ”) 
(70 7105 FORMATCIX, 15, 1X.F13.5,1%,F13.5) 
C871 «7110 FORMAT(’ VEIGHT VECTOR FOLLOWS: “/ 
| * INDEX“, 8X, WEIGHT’) 
(B72 «7115 FORMATCIX, 15, 1X.F13.5) 
B73 7128 FORMAT’ UNSICDESSFUL ADJUSTMENT WITHIN CRITERIA’. 
1 * PAD LIMIT CH ITERATIONS’) 
74 7128 FORMATC’ ADJUSTPENT SUCCESSFULLY PERFORMED WITHIN CRITERIA’) 
(75 «7148 FORMAT(’ ADJUSTMENT TERMINATED QUE TO COMMERGENDE ’, 


i “TG MINIMUM POSSIBLE ERROR’) 
EX 


OR RE ST 


om te 


4 = 
Qoeg 2 Grau Samed 
2 r) ® ' ® » 


eae 


Subroutine PWRITE 


HOIA-E «= TLE @-RY-61 83:52:43 
SUBROUTINE PYRITE 


C PURPOSE 15 10 WRITE THE RESAT FILE 


_ FORTH IY 
de at 
= ORD 
j 
<- 
| j 
+: OS 
i. ss 
; | 
, 
.. Oa 
fg? 
eee 
rece 
reG 
ml! 2 
watz 
(at? 


COPNON /RELOCK’ TRAHS, TSPPOR. OLIMIT. 
DLHILD.BEIGHT. 
POONPE .CODIPE . GAIN: 
POJINT, ERPAD! EPPLIM. 
> MON. WHIGH 
COPPIIN “IBLON/ TERA. IFRAT2. IFPPTS. 
TOLFID. NUPPAT . PTMSIZ. 
PROD] -NARAIT . IRDOHE: 
HATHT. TSS LENGTH, 
TRIG! . IBIG2. IBIGS 
DIMENSION TRANS< 2008) . TSPPNEC 200A) . 
DLIMIT( 488)  DUHILD( 402). 
WEIGHT (408), 


eine awn 


SPAN 


TFRET 1 (28), TFRAT2¢28), TFRMTS< 2a) 
IKSITMSIZ 
I=1 
(1 28 KEE TK 
TOET1+1TMS1Z-1 
WRITE(2, JFRMTS) (MEIGHT (1), J=11.12) 
T=1[+1TMS1Z 
CONTINUE 
RETURN 
END 


PRE i 


PSOIPE( 28) . CSOURE( 28) ,GRINC 28), IOLFHDC 28). 


! 
! 
| 


Table 9 


WTMAT. FOR 


Stochastic Adjustment for the "Map" System 


Purpose: The cornputer takes the transition matrix and 
provides a score matrix. 


Computer Output User Input 


e@ Run WTMAT 
<CR> 


Pause ~-- Insert data disc in unit 


one (DK:) and hit return <CR) 


*OK:TTEST1.DAT 
<CR> 
Name the result file 


*DK:TTEST2.DAT 


<CR) 


Enter the number of transition 
matrices (1) 


Enter size of row or column 
for matrices (I) 


I = integer 


erent en et te <r Ne ene OF 2 aT IRS I a vere 


Table 9 


WTMAT . FOR (Continued) 


Stochastic Adjustment for the "Map" System 


Computer Output 


Enter initial value of entries in 
weight matrix (RP 


User Input 


Enter allowable error in 
limit computation (R) 


Enter allowable error in score oe 
deviation (R) KOR) 
Enter maximum interations 200 
for limit computation (I) 
baceaiase as Ji SOR, 
Enter maximum interations 200 
for weight adjustment (1) 
<CR 
Enter format for transition 3 
5F3. 
matrices read (20A2) Sree 
Prooes oo S$ERD 
2R = real number, You must use a decimal point 


whenever the (R) appears. 


SWwhenever formats are indicated, parenthesis must be used. 


63 


Table 9 


WTMAT. FOR (Concluded) 


Stochastic Adjustment for the "Map" System 


wm. a = 


Computer Output User Input 


Enter format for score read (2F5.1) 
: (20A2) 
ae ee Tee ve = <CR> _ 


Enter format for weight output 


: (2042) (SF5.0) 


SoRORD tee, 


Enter lower limit for weight 
matrix cell (R) 


s 


Enter high limit for weight 
matrix cell (R) 


a Enter interation print switch 
for adjustment (1) 1 


tld Re came all a 


Example Output WTMAT 


Score Matrix 


i. _ FORTRAN IV HOIA-1 FRI @1-PRY-8} 2:28:88 PE 


Program Listing MLREG1 


C PROGRAM “MLREG! FOR’ 
C PURPOSE 1S TO PERFORM A MULTIPLE LINERP 
C REGRESSION AHALYSIS UTILIZING 
C SUBROUTINES FROM THE IRM SS PACKAGE 
DIPENSION X(68Q0) . XBARC 28). STD(28) .SEPC4AR) . 
1 RC218) SCPC 28) RECIE1) -RVCES). TFORMT 28). 
2 ROQEF C19), SDREC19), TUAL (19) ARSC cr, 
3 WORK! (28), WORKZ(28) » TSEC ZB) TMORE 19), NPD 1 
DATA X/6RRR+R. B. KBAR/2008. Be STD iE A 
i SCP 40GB. OY TRORNT ee’ FR ZIRHB, A SEP 2G. 
2 RABE 98 RY 1S RODE 1H, 
3 DRO 1SHEA, TURAL 1968, A DO, 
4 WORE 208, BY APR Ret, ISLE EG, 
i) FREAD. MNP AT, TER 2 1 ao 
PAUSE “IREERT (CATA DISC IM UNIT 8E CORE: af HIT PETES 
WRITEC7, 7118) 
CALL ROSIGHZ.' AE ROD CHS 
WRITE? 7115) 
READS. S116) 05 
WRITE(?7,7328) 
RERDCS. S116) NOPE 
WRITES, 7138) 
READS. S11@) NIURR 
WRITEC7,714@) 
READS. S128) CISAVECT) . T=1 HIURE: 
WRITEC7, 7158) 
READS. S11@) IDEPHT 
WRITEC7. 7164) 
READCS,511@) JUNIT 
WRITEC7, 7178) 
READS. SISA) CIFORNTCT) 324 
ISTOPSHOBS* (HURR- 1 +1 
DO 108 J=1. Noes 
READCZ, TFORMT CAT 3d, ISTCP NOES) 
IRERD=1READ+1 
ISTOP=1ST(P+1 
108 CONTIME 
CALL CORRCHOES, NAR.) BAR, STD SOP ROSE POH! LZ: 
CALL ORDERCHURR. R. IDEPT  MIVAR, ISAVE Re. 
CALL MIRUCRXNIURR, DETRIT » THNRKL TPES 
CALL PULTRCHOES NIVAR, CRAP: STD BSCP REP SALE. 
i RONEF » SORC. TURAL ARE. 
CALL ORDERCNURER. ICERYT HIUAR. ISAUE RABY" 
CALL LWRITECHOES. HUPE IVA ABRE. STD. SCP SCP RY PY TSE RCCEP 
SORT, TUAL. ARS. UNIT. IGEFIN 
STOP “EMD OF PROGEH MLREG!.FOR 
1088 FORMAT (32F8. 3) 
S11@ FORMAT (17) 
5128 FORMAT (1914) 
5138 FORMAT (2042) 
7140 FORMAT(’ NAME THE THRUT FILE cASSIaH) 
TINS FORMAT C’ ENTER HUMBER CF CRSERUATICS Gis > 
7128 FORMAT ENTER HUMBER CO UAPTRELES cis 
7138 FORMAT (® ENTER HUMBER CF OUADRHEOEIT CAR RELES C1)”. 66 


{ 
Program Listing MLREG1 (Concluded) 
| 
; FORTRAN IV HOIA-] = FRI SI-A-2) G2: PRE (22 
“WAT «7140 FORMATC® ENTER JEDICIES CF IMGERERCEHT URPIFELES 191499 
0842 «7158 FORMATC” ENTER INDE: CF DEPERDENSS UA RELE CID? 
_ OBS. 7168 FORMATO” ENTER THE CUTPUT UNIT HREEES .cLPs cSTTi> 
= B44 «7170 FORMATS’ ENTER THE IMP FORDAT SSEC FILA) C422} 
j 4S f ’ 
i e 
t 
' 


nee eT oe erento: 
a 


s 
’ 


ee 


ow 


ae : 4 ? 
eameng Date ei 
s 4’ « 
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Subroutine CORR 


PORTREN IV HBIFE} FRI @IHN-31 Gated PRE Gt 


SUBROUTINE CORR CHI 2. 5268 STDRE PED. 1) 
DIVEHGION KUL SRR TOU) GRP DR bb TaD 
00 168 J=1.¥ 
BCJ)=8.6 
108 1(J)=8.6 
Recent) Z 
DO 1@2 I-16 
102 R(1=8.6 


2922989228 
rs) -= 


. 
i 
—_ 


17 WHEEL) 
WARK DTC 

168 TH=TK1)-FH 
00 115 I=, 


DO 116 J=1.N 
Lau 
DID=KLI-TED 
NG BONDI: 
DO 115 Jett 
OO $15 KaI-J 
JK=JK+1 
1S ROKK MDC DE) 
Ike 
00 214 J=1.t 
HBPRCD SHEE CTF 
00 21G K-13 
K=IK+] 
218 ROKR -BODABE FH 


Z2B STOCI=SORT ORES RC IED? 9 
DO 238 J=1.H 
DO 238 K=J/8 
KCK AD /Z 
Loe I-14 
RACLERCI) 
LaecH-1)43 
PACLISRCK) 
IFCSTDCI)#5TDCK)) 225, 222, 2 
222 RK)=B.B 
@ 10 28 
TB ROMER STO I ASTOR) 
238 CONTINE 
FHESORT CFF 1.6) 
DO 246 J=3,1 
248 STOXCD=STDCI AN 


a | 


Subroutine CORR (Concluded) 


d. —pORTRNTY | HBIT FRI G1-tHH-8 wns) PRE 082 
j oss) Lats 
b. 6  BO BDL) 
on? RETURH 
( RE B® 
| 
| i 
cj 
j 
j 


+e EET. $8 Ue 


Subroutine ORDER 


FORTRAN IU Halfet FRI GI4R-e! 62520059 PRE @i 


QBROUT IIE CPDERCN RO HDEP. K IAL BOR: 
DIMENSION ROL), ISMECLS POU Re? 
EO 
00 138 JeL.K 
L2=ISMECS) 
TFQMDEPL2) 122.122,125 
122 LAWDEP+ (L2H 2422/72 
@ 10 13 
123 L=L2+CHDEPHDEPADEP }7 
1B RD=RL) 
00 138 I=1.K 
LISISAUEST) 
IFCLI-“L2) 127, $28. 128 
427 LALI4CL 22-4122 
@ 10 123 
$278 LA2eCL PL I-L 1/2 
129 WEHH! 
124 RACRDERL) 
ISE CKD DEP 
RETUPH ; 
a, 
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I Subroutine MINV 


¥e 


2 FORTRAN TU MBIA FRI GEHIR-et GitzZ1212 PE OE 
= Bt = ABROUTIFE IMU DL 
{ez iret nC) Le ters 

eS .E 
eH 
} eh 00 88 K=1.N 
MRK 
(7 LK 
BBR KKH 
09H 
ere BIGAFAKKIO 
11 OH 
tz ICH 
1s Om ISH 
1g AIZH 
OBIS 10 TFLABSCRIGAY-PASKACIDD? 15,2028 
e816 15 BIGRACI 
@17 Wl 
eis KK= 
019-28 CONTIN 
we JALAK) 
210 IRD TE 
_. 2 BAH 
23 008 I=L 
- 24 KSAT : 
5 HU HKD) 
me TNSRI-KH 
wze AKI =AKIDD 


we 863 ACIDE 
os 8239 3B ik 

CS IF) dh.45.28 
~ C31 38 JPeCI-1) 
-@2 00 40 JF1-h 


03 KAR) 
084 jI=JP+] 
835 HOLD= =F) 
=) BS AKA . 
37 BW ACDHLD 
3- @BBB «45 IFCBIGA) 48.46, 48 
i. @9 «Dee 
040 PETUIRN 
- RL BOS Itt 
{ —etz_IFCI-K) 50,55.58 
084350 IKKE! 
B44 ACIK)=ACIK)-<-BIGAD 
| 0S «SS CONTINE 
6046 00 65 I=1,H 
07 1KK+1 
| ese Weis 
649 00 65 JeLN 
- 058 LelJeH 
0851 IFCI-K) 66,65, 68 
{ O52 IFUK) 65.02 M 
“@S 62 KJs1J-14K 
54 ACIDEAC TOAD HLS) 


yy I Subroutine MINV (Concluded) 
i. FORTRAN IU HBIFRL «FRA GIPRN-B1 @28Z1812 PHE 22 
j 55 65 CONTINE 
- KKH 
0057 00 75 Jett! 
; 5B KIEKIHN 
{ee IFCI-K) 70.75.78 
7 EB TB AK=AKD-2108 
_ RBs 7S CONTI 
j 0862 D=DYBIGA 
-. ORS PKK) =1.0- RIGA 
0064 9-88 CONTIN 
i 065 KH 
| - O66 108 KK-1 
067 IK) 158.156, 185 
; eee 16 11) 
{ ee TCA) 1 120. 1 
C878 108 JONeCK-1) 
a7 aed -1) 
72 DG 116 Jett! 
: 873 KeJ0+) 
a74 HOLD=AC JI) 
Sc) NeIR+] 
~. | are AIK)=-ALII 
P 77 | 11G ACID HD 
-. OR 18 JAD 
ars TIA 1. as 25 
~~ 1S KEK 
| 00 134 I=1.K 
wee |. <44 KIKI] 
> HOLD=ACET) 
84 JKI-*K+) 
= eS AKD=AID 
$. 6886 138 AUD 
087 TO 1¢2 
(ORS «= «158 CONTIIFE 
i pees PETE! 
en END 


Gaatiae 
¢ y 


rs 


Gemeq iomsonen @mamany 
- . . ; 


' 


Seg eyes Rie ggg ce mee 


Pi ere are, ae rite 
3 ‘ e ‘ 
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Subroutine MULTR 
WOiAe 1 | FP] Gifs 81 Ges cit ci PAE 1 


SUBROLITINE PULTRER.E ERS. STL LSE B56. DS 


OIVENG ION SBARC1) STOCED GLE Rady abd 
i SAE BSE Td tc! 


EK+1 
00 108 FI.K 
108 BC J)=8.8 
DO 118 J=1.K 
Liske<J-1) 
00 11 I=1.4 
Lal+! 
118 BOD=BC DHCD PL? 
60-8.8 
LI=ISAUE CHD 
00 128 I=1.K 
REED 1D) 
L=ISHEC]) 
BOD BCD (STIMEDSTIOCLY) 
128 BOFBOHB( DEAR CL? 
BORKBAP (LI): 
SOPR=PIMANL I: 
PPOPT PES PT} 
SSOR=DEL I-58 
FRSEK-1 
SY=SSOR- FI: 

00 138 J=1.k 
Liskec I-14) 
=SAECT) 

1F BODO ARERR AML 5 
138 HYBIRT 
135 SY=SOPT CPR SI 
FKK 
SSPRMESSAP FT 
SOORMESSOR. FH 
F=SSARI- SSDI, 
fG61)=80 
FASC 22% 
FASC 325 
FAS 64) =S08P. 
FRS(S) =FK 
FS 6)=S5FH4; 
FHEC7)=5DR 
ARS(8) =F 
FNS(9)=S2DPH 
PRS (18)=F 
RETURN 
£0 


fd iii ae, 


2" Subroutine LWRITE 


Lm) 
FORTRAN 1 HAP FRE G1teH~21 Ges zisa7 PRE oS, 
i eee pees LEPITES TPERD. HENLE. BEE. sm. SF, OSCF. aa AMEE FS na 
eS POOEF OPC. TL, Ta tO 
ene ‘HENSIOH YOR). STEEL SAP LS, DREADED) 
1 SPEC) PER DPT TL RRL 
£0 NURRESNIURR+ 1 
one4 WRITECIUNIT 6110 
0005 WITECIUNT 6120) IRE 
, eRe WRITECIUNIT 6128) FRE 
ne a WRITECIUNIT 140) IESE 
is) WRITECIUNIT.G156) CISH#E D3 121 HSE: 
Ci) WRITECIUNIT 6160) [SA Sc eest : 
wie WITECIUNIT. 617%) 
eet! 00 108 I=1,HIUeP 
12 JEISMECD 
e213 WITECIUNIT.618 IstoecD Stic: 
0814 108 CONTIME 
015 PRITECIUNIT 6190) YROPCTHOEP' sTbe Tee” 
ate WITECIUNIT.6191) 
7 00 116 I=1,1uFe 
12 00 165 J=1. Hert 
eng HEISEI 
>. O20 JETS) 
0021 IFCULLE. LL) Gi TH 108 
9823 KeCLi-1 Hee 
Bt @n24 QRITECIUNIT. 61923 TT. SOR 
0225 «105 CONTIME 
~ 082% 118 CONTIME 
+ 27 WITECKAIT 619 
4 RAZR D0 128 T=. Hike 
eas DO115 Jey Mie 
> RB IFQ.LE.1) G75 115 
we» BZ K=CJ—1 104} 
3 WRITECIUNIT 6194) ISHED) ISHECI PS 
7 0054115 CONTINE 
wa, 8S 128 CONTIN 
ORE WRITECIUNIT 6198) ISHESHEL) 
~. CW 00 125 I=1,Hiee 
; BB WRITECIUNIT. 619) ISAECT) Fact 
“- @39 125 CONTIME 
ete WRITEC UNIT, 6197) 
7 gaa DO 132 121, NURR: 
$. (O42 J=ISAUEC1) 
ro 43 IFCDSCP( 3) GT .999999999,8) GO) Th 123 
sm - O45 WITECIUNIT.6198) J.0609C 1) 
i 146 @ 10 1% 
. 0047 «123 CONTIMLE 
_ 0048 WRITECIUEIT.6193) J.DECPCD) 
{ 6849138 CONTINE 
-- 9B WRITECIUNIT, 6283) 
> 51 00 208 I=L, HiUde 
3° 2 WRITECIUIT 6210 TSAECD) POLED) Sher TRL) 
7 053-28 CONTITLE 
054 WRITECIUEIT.6220) FEC LA 


- 885 WITECTUNIT 6238) FED) 


=, -* Subroutine LWRITE (Continued) 


~ FORTRAN 1U HBIA-] FRY GIS Gzezisd? PRE G32 


wre ROSFEECZ)842 
87 WRITECIURIT.6251) RZ 

WB WRITECIUNIT.6248) FEES) 

0259 WRITECIUNT 6255) 
en WRITECIUNIT 6268) AREA) ARGC) ARECG) AE C1) 
0261 WRITECIUNIT 6270) FEC?) ACR) ECS 
0862 RETURN 
| | 5 B63 GIG FORMAT’ “7° 7° MULTIPLE LIFERR PECEESSION) 


Bw ae. t 


" @B64 6128 FORMAT’ HUMBER OF CRSEMATICHS = 017) 

L @865 6138 FORMAT(’ HUBER CF URPIPELES = “0173 

i OB65 «6140 FORMAT.’ HUMBER CE IDEPEIIEHT WHR TPELES = 1 

: : 0067 oar ale WDICES CF INDEPEMIEIS LMP TFELES: 

* 4914) 
: @B68 6168 FORMAT(’ INDE: CF OEPEMDERT UHETPELES oD 
: ORES GI7A FORMAT’ *-° HEARS FAD STADPPD CEU IGE GF TERENCE 
“URRIRELESS 
* INDEX’ Qe. (HEAR au, “SDELE - 


ee OR oe mH Soe, 
reyes sans 


by ot 


wey 


0678 6188 FORMATCIX ISIN FEL. S LR FIL DD 
@B71 «6198 FORMATC “7° EAN OF DEPENDENT URPIRELE = OFil.2 
4 Pes ! * STANDARD CCUTATICN OF DEPERTIS USGrELE = Fil 
\ OB72 «6191 FORMATS: 7" SUN OF CPOSS PROT: 
! “INDEN  2 INDER ere 


? 


$3 


Fy oe ee 


. 
rast 


; 2. ee 
} = PBFS 6192 FORMAT CAN, 15.20. 15-18.F11.2) 

8749S FORMATO’ 77° INDEPENDENT WPTEELE COFHELAT OH 
17 TE 2 TIRE eo 


| ee 
3 ¢ ‘y 
we OBS 6194 FORMATCIX. 15, 28.15, 1279.3) 
; @876 6195 FORMAT(’ “7° CORRELATIONS WITH DEPENDCHT “Se iRELES 17 
- ae 1 * TNDEX’ 6%. “COPR’ 
: et 2 imei. ieee 
3 wey 
= 0877 = 6196 FORMATCIX, 15. 1%,F9.3) 
; B78 «6197 FORMATC’ “7° SM CF SURES: 
pa { * INDEX’ 9%, “SUMS 7 
2 - a aes 
: 3 77) 
0879 «6198 FORMAT(1X, 15, 1X,F13.3) 
oe 0863 «6199 FORMATC1X, 15, 1X,F13.0) 
OBB 


6280 FORMAT’ ’7° REGRESSION COEFFICIENTS: 
: * INDEX’ NOR, COEF’ 10H, “SDEW’ 5, “T-URLUE 


1 
2 Fn 1 I eS 
3 pe 
j OBB2 «6716 FORMATCAX, 15, 1X,F13. 7, 12,F13.7, (X-Fi. 3) 
‘ - B83 «6720 FORMAT’ ’/’ INTERCEPT = “.FI1.3) 
* @BB4 6238 FORMAT(’ ’/’ MULTIPLE CORRELATION CORFFICIONT = ‘.Fe.2) 
{ B85 6731 FORMAT(’ ’7’ EXPLAINED URRIARCE = “.F6.3: 
i _ OBR 6248 FORMATC 77° STANDARD EPROP CF ESTINATE = -FII.3) 9 75 
B87 «620 FOPRAT(’ “7 FRALYSIS (CF UPPIARTE TRELE ~ 
| * SPURTE’ 2X, SUINGIY ES. DE a. TEORSO RE. F-Pe 


Subroutine LWRITE (Concluded) 


La 


- PORTRAR 10 HBIP-1 FRI GI-PRN-BN 2521247 PRE G53 


2 a SR a es 
@. 3 4’ fy 
(083 6268 FORMATO: REG’ -AX.FIZ. IPR PE SS PLZ 2) 
0289 6278 FORMAT’ RES”. 4K.F12. 3.17. FG.G FITS 
0898 8D 


LTR PRS Cn A er BAL, EEA AS Geren = 


Ty 


Table 10 


MLREG1.FOR 


"ore Multiple Linear Regression 


Purpose: Performs a multiple linear regression, 
Explanation: 
Enter # of observations = total # of observations 


Enter # of variables = total # of independent and dependent 
variables in your file 


ee a. 


Enter # of independent variables = total # you are going to 
use for this regression 
Indices of independent variables = the number's that correspond 


to the independent variables in your file 


Computer Output User Input 


e Run MLREG1 


Pause -- Insert data disc tn 
: unit one (OK:) and hit return 


<CR> 


Name the input file (Assign) *DK:F 20935. DAT 
<CR> 


Enter number of observations(1) 180 


<CR> 


Enter number of variables (1) 16 


<CR> 


77 


Se ee a Ae Cee Se 


. Table 10 , 


MLREG1.FOR (Concluded) 


Multiple Linear Regression 


Computer Output User Input 


Enter number of independent 
variables (I) 


Enter indices of independent 
variables (1914) 


Enter. index of dependent (I) 


Enter the output unit number 
(6=LP, 7= TT) 

<CR> 
urmter the input format (16F 20. 10) 
specification (20A2) 

¢CR> 


1 . 
Whenever formats are indicated, narentheses must be used. 


Example Output 

RA TIPLP LINFAR REGRES ISM 
WIMFFR OF OBSERUATIONS = 174 
WIMFFR (F VARIABLES - 16 
MUMGER OF “JHDEPEICA HT (MRIABLES = > 
HOICES OF JNDEPENDENT VARIABLES: 

a 10 
INDEY, (CF DEPENDENT UAPTHELES 36 


VERE AAD STRDRRD DEVIATIONS OF SEP ED ENT Une ABLES: 
THEY. VFA “OEW 


8 0.297 6.292 
16 6.472 0.425 


PEAR CF DEPENDENT (HRIAFLE = 6. 4s 
STRADARD DEVIATION CF MEPENDENT URRIAELE = Q. 38) 


et CF CROCS PROOUCTS: 


THDEZ.  THOE:: IWP 
% in 12.429 

a Ir. 1S, oe. 

sa fr. 12. Gra 


THDEPFHOFT VARIABLE (OPRELATICNS: 
THORS. TES: CORR 


re | SS 


(TORE ATIONS. UITH DEPENDENT VARIFELF: 16 


IEY (rep 
2 4, Oe. 
1@ Gerd 
QM OF CURES: 
JHDE:: AI: 
2 it enn 
1a TIS. 
16. Fie ( 


REGRESSION COEFF IC TFNTS: 
INDEX COEF SEU TALE 


_—— —_ —— 


g a, TMX. B16 R284 i] | othe 
i a cs P| 8, Ode) 32? tedle 


INTERCEPT = 8. RAS. 


MAR TIPLE CORRELATION CCEFFICIENT = Q.250 


MLREG1 


Example Output MLREG1 (Concluded) 
. > 
FYPL AIMED URRIAE = a7 


STANDARD ERROR CF ESTIMATE = 244 


ANALYSIS OF URRIFOICE TRELE 
SURE QCA ME rec Far iv 


ws eee ———_ 


REG 17.569 : 3. Mey. 
7.5)! 8.04552 


| 
7 
L. 


eer amen mem me 


= oe, 


I Program Listing HADMRD 
+ FORTRHA JU Helf-i = FRI. @-HAY-21 08238: 48 PRE ei 
‘i © PROGRAM. “HADPRD. FOR’ 


(: CREATED FOR THE HINA: 9-JAN-86 
f PURPOSE 1S 10 PERFORM A FAST HADAMARD TRANSFORM 


c C FOR A URRIABLE NIMBER CF BONLEAN TIME SEQUENCE RECORDS 
c ( CONTROL INPUT DEFINES THE SIZE OF THE HADAMRPD MATRIX, 
( NUMBER CF BITS IN BOCLEAN TIME SEQUENCE. NUMBER. CF 
i ( SHIFTS TQ PERFORM CH BTS, AND THE LENGTH OF SHIFT 
i C THIS PROGRAM 1S A PART OF THE ARSOLLTE COePLTATICH 
=: ( PROCEDURE FOR THE “MRP. SYSTEM 
0641 COA /ABLOLK/ 18. THDDIN, IPCAD 
f (enerz DINENSICH TETSC1624), TRC128) STATC1 28,29 
ES. PRISE “INSERT DATA OISC IN UNIT CHE CDRS) APD HIT RETUPH’ 
CT MPUN CETL UPR IAELES 
ed WRITE? 7108) 
jae CALL AEST 1. POO HED) 
(hts. WRITEC7 TH 1G? 
a7 CALL ASSIGH 2." .-1. HE NOD 
abet WRITEC? 7120) 
a 5) PERCE SiG) TENS 
ese WRITE? 7138) 
GN PEAS. SIAR) LERETS 
NZ UPITEC? 7140) 
AIS. PEAS. S109) HSHIET 
etd IPITEC7 7158) 
cS bal PEAS. S108) LENSHE 
a ( COMPUTE OTHER COETROL UPIRELES 
PALE THOME TPO 
] i PSIZESFLCAT<THDOIM) 
(ye NOME L=HCHIET +4 
my 4 FESIZE=FLOGT “HOME! ) 
~ Met IPERD= 
a ( PERE (HE BTS PECL 
Met fete COAT 
qe PEA L J108-END=7AQ) CIBTSCT), I=} LEMETS? 
NM (7, IPEPD=IREAD+1 
(tg IFCLENBTS, CE. 1024) G10 128 
(fe, ISTART=LERBTS+1 
1 7 M0 11@ J=ISTAPT, 1024 
(are, IBTSCI)=@ 
(K79) O18 COATINLE 
we ER CONTIMLE 
j ( COUNT HUMBER OF ONES AND ZEROS IN BOOLEAN TIME SEQUENCE 
° Me ItE=@ 
aan? IZEROR8 
[ a, 10 138 }=1 LENBTS 
0 4 JFCIBTS(1).£0.1) JOMES1OHE+1 
(0, IFCIGTS(1).EQ.@) IZERO=IZERO+1 
- C% 120 COATIME 
[ ( INTIALIZE POPULATION TOTALS 
(a PHTOT=0.6 
1 aes PHTOT2=0.0 
(a4) PUTOT=8.8 
gy PUTCT2=8.8 


( INTIALIZE STATISTIC MATRIX 


81 


ao 
I 


— sth 


Program Listing HADMRD (Continued) 


| KETEAA I! «= HOIARL «FRI GB-MPY-1 @Bz 38248 PAE 02 
(43 Xt 200 321, IHDDIM 
T waa STAT CL 1)=8.€ 
me gh STATI 226.8 
se 2001 CONTINE 
: C SHIFT BTS TC WORKING VECTOR 
7 DO S60 T=1 NSH 
ie: INDEN=@ 
_ ge JSTABT=(1-1 #LENSHE +1 
Ax JSTOP=ISTART+IHODIMHL 
7 Gel 101 308 JRISTART. JSTOP 
wee THDEX=THOEX+1 
SS IRCINDEX)= BTS) 
L. (itd TFCIRCINDE).EG.@) IBCINDEX)=(-1) 
OK Et CONTIME 
¢ PEPFORM FAST HADRMRRD TRANSFORM 
enr CALL HADES 
C GATHER. STATISTICAL DATA 
(es 00 402 J=1, JHDDIM 
wea STATI, 1)=STATCI, 1 )4FLOAT CIBC) 
(ees STAT (J. 2)=STAT(.1.2)4FLOAT (IBC J #282) 
G1 ae CONTINLE 
Me? SE COTINLE 
COLTPUT THE SOLUTION FOR THE CURRENT BOOLEAN TIME SECUENCE 
(ne3 WRITEC7,7168) TREAD, LEHBTS, 10HE. IZEP(. THDOIM, HOHIFT. LENSHE SHE 1 
4 WRITE 2,210) TREAD. LEMBTS, IONE. ZERO, THOOIM, HSHIFT LEHSHE. HEHE 1 
ae 10 628 I=1-JHOOIM 
* Mttes PEPHESTAT (I, 1)RSSIZE 
ee? IFCHGHFI.LE.2) SUPP=0.8 
(ne IFCHSHF1.GT.2) SURR=(STATCL, 2)-(STATCT S9STATCL, 1)/RSSIZE)) 
a 1 ARSSIZE-1 
er! ERITEC7, 7178) 1,STATCE, 12) STRTCL, 2+ SERN, SURP 
- rz UPITECZ,2110) 1, SMEAN, SUR 
-. OS PHTOT=PHTOT+SHER 
mars PHTOT2=PHTOT2+SE RH SPERH 
-- #5 PUTOT=PUTOT+SUF 
A PUTOT 2=PUTOT2+SURRASURR 
=" 7 6H COMTIME 
Pe: PHRERHEPITOT PS1ZE 
en te) PAAR (PHOT 2-<PHTOT#PMTOT PSIZE) )/PSIZE 
or 6M PUERPERPUTOT PSIZE 
(eee PAURR= (PUTOT2—(PUTOTSPUTOT PSIZE) ) PSIZE 
~ Me? UPITE<7,7188) PPMERH, PHURR, PUMERH, PURE 
7 ees WRITE(2,2120) PHMERH, PHURR, PUMERH, PUJRR 
( PEAR NEXT BONLEAN TINE SEQUENCE 
22) Gi T0 1% 
i C ENE CE PRORRI 
ees IL: CONTIMLE 
C2 «= STOP “END OF PROGRAM HRDPPD. FOR’ 
| Ce? 160 FORMAT C1281) 
1G FORMAT (BIS) 
4 4 NG FCPRRT 13.269, 3) 
qe AS FORT CAFI3.3) . ' 
| EAS She: FCOMATCI7) 
Cosy 7100) FOPRATC’ HARE THE INPUT FILE CRSSI@N’ 7) = 
[. 


THA I'S HEIL =—- FRI @OHIRY-OE G38: 48 PRE 3 


a Program Listing HADMRD (Concluded) 
© | 


re 093 7H1@ FORMAT<’ HARE THE RESULT FILE <ASsI@!’”) ! 
Bees 7I2G FORMAT <’ WHAT PCKER OF 2 EQUALS HRDPRRPD SIZE (1)’) | 
5 7138 FORMAT’ WHAT IS LENGTH CF INPUT BONLERN TIME SECUEME (1)”) 7 

7 (9S 7146 FORRT(’ HOS PRAY SHIFTS ARE TO RE PERFORMED CH SPPPLE (1° 

3.0087 TISE FORMATC’ WHAT. IS THE LENGTH CF A SHIFT <1)’ 

RSE 7168 FORMAT(’ FAST HADAMARD TRANSFORM’ / 

* INPUT BTS = "17 

* LENGTH CF ATS =. 577 i 

* WIRRER CF ORES TH BTS = 7,17¢ 

* WEEP CE ZEROS IH BTS = 7.177 : 

* HADROPRD SIZE = 4177 

“ EPREP CF SHIFTS =7, 17° 

*” LENGTH C& SHIFT = %,177 | 
| 
i 


Came. 5 


* SHPLE SIZE = ed 


° HRORRRDY 7K, SPR SK, “SLPS 7X, “PERN 1, 
UPIRLE’ 


f rhe i eK aK, ee, 

ns 
O99 IPR FOPMATC’ °01%,17. 2012, F 18,8), 2(2K,F9.3)) 
CIE «71S FCRNATC’ POPULATION STATISTICS’/ 

\ * WERE! CF FERNS = OFIRY 

Zé * URPIARDE (F MEANS = * FIL 

3 * MEPS OF URPIFATES = “FILLY 

4 ’ URPTANCE OF URRIANIES = “,F13,3) 


7 op Oo DONEMEUNNNO 


FOP TREE 1 HBIA-] «FRI. GBF! 08:59:69 


SUBROUTINE HADMRD 


C PURPOSE 1S TO PERFORM A FAST HADAMARD TRANSFORM 


C (tt At! INPUT BOOLEAN TIME SEQUENCE 
COMMON ABLOCK/ IB, JSTOP. ISTOP 


z 


—~ 


$4 


4 4 Fb! 


| 


v 


Gtith tomnand pony | 
« e a ‘ e ‘ 


‘Qumed 


v 


ri 


DIMENSICH J(128), 1BC128) 
00 KE =f, 1ST0P 
INDEX=1 
JINC=2+41 
KSTOP=204(1-1) 
DS 388 J=1, STOP: SINC 
DO 108 K=1,KST(e 
LOC1=J+K-1 
LOC2=LOCI+KSTOP 
JACINDEX)=1B¢LOC! +IB(LOC2) 
THDEX=TNDEX+1 

100: OOUTINLE 
DO 20% K=L RST 
LOCIEI4+K-1 
LOC2=LOCI+KSTOP 
TACIBDEX=1B(L0C1 )-IB(LOE2) 
INDEX=INDEX+1 

222 COUTINLE 

RE COMTINLE 
(M1 468 J=L, ISTO 
IECDSIACI) 

S00 COTIMLE 

S43 CONTINUE 
RETURN 
EN 


Subroutine HADMRD 


Program Listing RLOGIC 


| Pree 4d 


“FORTRAN! HAIA-} FR] GR+HAY-R1 688 29:49 PRE 68! 
C PROGRAM “RLOGIC. FOR” 

€ CREATED FOR THE HIJA SYSTEM? 8-JAN-88 

( PURPOSE: ACCEPT A FILE OF BOOLEAN VARIABLES 

C AAD A LOGICAL GOAL VARIABLE. COUNT THE NUMBER 

(CF OCCURENCES FOR ALL POSSIBLE COMBINATIONS 


( (CF THE MINTERNS. COMPUTE THE PROBABILITIES 


(. 
¢ 
¢ 
¢ 
if 
{. 
( 


¢ 


FOR THE EVENTS OF ($) MINTERM TRUE 
<2) GOAL TRUE GIVEN MINTERM TRUE 
C3) GOAL FALSE GIVEN MINTERM TRIE. 
TEST FC@ PURE LOGIC FUNCTION. WRITE CUTPUT 
YO TERMINAL. WRITE CLTPUT FILE OF REDLTS. 
THIS. PROGRAM 1S A PART OF THE RELATIVE 
TCHIC PROCEDURE FOR THE “MAP SYSTEM, 
COON -AACOTKT PRO. ROUNT .KUALUE. PRRADT. 
TRGANP LOGIC, HURR: IPURE- LENGTH NSIZE 

DIMENSION PROB 324, 2) .KOUNT(1624,2), KUALUEC 1024), IBSAHP C1) 
PAUSE “INSERT JMPUT DATA DISC IH UNIT CHE (DK:) AND HIT RETURN’ 
HACE TAPUT FILE 

WITEC7 7108) 

CALL RESIGHEL.° 74-1. °R00" “HED 


J 


C HACE CUTPUT FILE 


WHITE? 711@) 
CALL AGSIGH(2,° “.-1, NEW, HOD) 


( ENTER HUMBER CF URRIABLES IN SAMPLE 


WRITES? 7128) 
PERDCS, S1@8) HUAR 


( FITTER THE COEFFICIENT ADJUSTMENT FACTOR 


tRITEC7 7138) 


ait PEADCS, 5110) PREADJ 
( INTIALIZE URRIBALES 
ante: IPURE=@ 
oie HEIZE=R 
mtd LENGTHS 2hHHURR 
eats 1X 106 1=1, 1024 
MSS YOUN, 1928 
15; HUNT CL, 258 
Coe KOUNT1,3)=8 
i (ang PROECT, 1)=0.8 
Pe aK PRIECI, 2958.8 
eres PREBCI,3)=8.8 
> an HURLUEC1)=-1 
{ (607. 1F(1.GT.18) GO 70 168 
ws IBSAPP(1)=8 
- WO 10 COMTINE 
i € EAD ONE RECORD AND COUNT 
> Me? 7 COUTIME 
L car NGIZESNSIZE+! 
ran CALL BTALLY 
#2 @ 10 28 
C COMPUTE PROBABIL TIES 
(57 28 CONTIME 85 
can CALL PLPROE 


€ TEST FOR PURITY CF LOGIC FUNCTION 


+ ee 
“ 


i Program Listing RLOGIC (Concluded) 


- SCR TRAN I HOIA-1 «FRI @B+HRY-81 68:29:49 PAE 82 


+ 004 CALL TSTPUR 
¢ MRITE 10 TERMINAL DEVICE 
ans CALL TWRITE 


ee eee 


2 C WRITE 10 FILE DEVICE 
MARS PAUSE “INGERT OUTPUT DATA DISC IN UNTT CHE <OK:) AND HIT RETLPH’ 
E27 CALL FURITE 


| Oe STOP “END OF PROGRAM RLOGIC. FOR’ 
39 1G FORMAT CHIT) 
40 (5108 FORMATCI7) 
© RAS | SIG FORMATCFI3. 7) 
* 47 7100 FORMAT(’ WARE THE INPUT FILE (ASSIGN) /) 
C4 | 711G FORMAT(’ NAME THE RESULT FILE <ASSIGN’/) 
M44 7120 FORMAT(’ ENTER NUMBER (F URRIABLES FOR PROBLEM <1)’) 
M45 IS) FORMAT(’ ENTER PROBABILITY FACTOR FOP *. 
) “COEFFICIENT ADRISTHENT RY’) 
ite. BD 


: Subroutine BTALLY 


{ Wet SLEROUTINE BTALLY 
: C PURPOSE 1S 10 COMPUTE LOCATION FOR ADDITION OF 
{ CUPPENT SAMPLE TO TOTAL AND THEN PERFORM ADDITION. 
ORE? COMMON AALOCKI/ PROB. KOUNT KURLUE, PREADJ, 


* FORTRAN IY! HOIA-1 = FRI GB+HAY-@1 08238204 PRE 981 


i ! IBSAPP LOGIC MURR, IPURE, LENGTH, NSIZE 
QE. DIMENSION PROB 1804, 3) ROUT 1@06,3), KUALUEC1€D4), TESPRPC1@) 
HDI 
] Geet D0 108 T=1,NURR 
* 20% = TFCIBSARPCI),E0,1) TNDEXEINDED200(1-1) 
PER 50 CONTINEE 
| | RG FOUN CINDER, 1)-AOURTTCIRDER, 1941 
LPN TFCLOGIC.E0.1) ROUNTCINDEX, 2)=KOUNT CIHDEX, 2)41 
@R12 ———_TFCLOGIC.NE. 1) KOUNT THER, 3)XOLNTCINDES: 3)41 
| m4 RETIN 
| eis 


i 
; 


Te ee  F 


Subroutine RLPROB 


FORT ROE 3! HOIA-1 FRI @B-+HRY-A1 Gas fe. PE OO 
att} SECUTNE PLACE 
ie C PUREE 1S 1G COMPUTE THE PROGARILITIES FOP LOGIC FLORTICE 
(enirz COA ARLOCKI/ PROB. KOUNT. KUALUE, PREADT. 
! TESAMP LOGIC. NUR, IPURE LENGTH .HSIZE 
. (G. DIFENSIOH PROG 1624, 3) .KOUNT( 1024. 3). KUALUEC 1824) . JECAMPC1G) 
a) PSIZE=FLOATCHSIZE) 
mr DG 1@6@ 151, LENGTH 
, (ete. TFKROUNTCL, 1).LE.@) Gi TO 168 
hd COURT 1=FLO@TCKOUNTCT. 1) 
(an COUNT FLOAT CKOUHT CT. 2)) 
Alt COUNT 2=FLOATCROLUT CL. 2)) 
mals PREC TD =COLMT 1 RG1ZE 
i PCR) 2008.01 2- UT 
(e117. PRET. =COUMT VCMT 
(Old tet CTI 
fais ETL 
ie. ED 


oe 


- emerge 


s- Subroutine TSTPUR 
=m 

FORTRAN T= HAIREL FRI @Q+MRV-@1 08238229 PROE. Oat 
= A QUPOTINE TSTPUR 


C PURPOSE IS 10 TEST FOR THE PURITY OF THE COMPUTED LOGIC FUNCTION 
RAR? COMIN ABLOTKI/ PROB. KOUNT KURLUE. PRERDJ, 


: ! TSSAMP LOGIC, MURR, IPURE, LENGTH, NSIZE 
i Ss DIMENSION PROB 1624.3) KOUNT< 1024, 3), KUALUE (124), IRSAMP(1@) 
ead 0 168 1=1. LENGTH 
GER TFCKOUNTCT. 1).LE.@) KUALUECID=@ 

* eR? IFCPROBCS,2).GE.PREAL) KURLUECD=! 

(nr TFCPROBCT, 3). GE.PREADI) KUALLE(1)=8 
al TFCKURLUE(]).LT.@) G1 TO 108 

sts IPURE=1PURE+! 

Ald 101 CONTIMLE 

i be RETURN 

CALS EN 


89 


inten eS ot i ee ee Meee deme ee 
Se kart 
: wd .- ai 

: ‘ a iat ta. wt ates male 


. 


= 
Subroutine TWRITE 
~ FORTRON = HBIPEL «FRI @Q+HAY-81 Ga: sae4! PRGE O81 
Caen! QUPOUTIKE TURITE 
C PURPOSE IS TO WRITE TO TERMINAL 
een cre /RLOCKI/ PROE. KOUNT KURLUE, PREADJ, 
TRSAMP, LOGIC, NUAR, IPURE, LENGTH, NSIZE 
GRC eee PROB 1224.3), KOUNT (1024, 3), KUALUE (1024). IEGAMPC 103 
eeeteg HPITEC7. 7108) HURR. LENGTH, PRBRD 
GCS VPITEC7,711@) NSIZE. IPURE 
GOK, ITEC, 7120) 
eA? DO 10 {=f LENGTH 
paee WPITEC7,7130) T.FOURTCL, 1) .PROBCL, 1) KURLUECT, 
{ KOUNT<1.2),PROBCL.2) KOUNTCL.3) PROBE, 3) 
Rs 168 COTIIE 
pete PETIEN! 
GCl! PGR FOONST.* RELATIVE BOOLEAN LOGIC SEARCH FOR A LOGICIAL ’, 
"PUMCT ION C& THE MINTER’ / 
2 * HPTEP CE WOPINES = 177 
3 SURMOCE CE NIETERNS = 177 
4 " OMCTNENT FACTOR = 70F9,5) 
PAZ 1G FOOMATC IPeER Ce come. = «017. 
j EIT = IPs 
7: PITT SHN®D ECU THE HUPRER CE MINTERHS ”, 
4 ee Fe MEE LOGIE) 
eet Yt DEOURT 5 “NPE 
| re 


OAl4 FILS FoereT ©: Ges piesa 1 eRe LE ea 
! my tt 


a8) fal ef 


90 


3 
E 
Le 
fe) 
> 
q 
Oy) 
1) 
ie} 
4 
O 


v 


tm 4 


Se ed el 


| a 


Subroutine FWRITE 


‘ 
: F rorernn 1U HMR «FRI OSHPN-BI 023053 PRE O8t 
# Cet SUBROUTINE FLITE 
| : [ C PURPOSE IS TO WRITE USTPUT FILE 
+ Liew CORRE /ALODE'S/ PROB, H'OL2ST  KURLUE. PREAD 
e : 1 TESATE, LOGIC. HURR IPURE, LENGTH, NGUZE 
a Pg es DIMENSION PROBC1224.3) KOU (126, 3), KUALUECL@24) » IBSRPPC 18) 
4 [_ wm 9004 WRITE(2,2108) HIRE. LENGTH, PPECDY, SIZE, PURE 
D0 108 1=1,LENSTH 
po WRITEC2,2110) LLRCETCI 1) PROBCE 1) FURLUECD, 
[ ' ROR (1, 2), PROBCI» 2) HORSTCL,3) PROBE 3) 
6087 108 CONTINLE 


HM CORI = 2108 PORPETC217,F13.7,2177 


‘ f OBI 2198 POPPAT(217,F8.3.13,37,F8,3.17,FS.2) 
it 81! BD 


Gece oP, —! 


é 
® 


hs 
PEER a8 oy ageneton- eae 


Guam 
et 


ee 


ae | 


; 8 
ALG 


GLOSSARY OF TERMS 


ie <CR> Carriage Return | 


eae wad 


LP Line Printer | 
aa | 
{ TTT Teletype 
| 
re) Integer 
3 
®) Real number. A decimal point must 
|. be used with real numbers. 


tng 


NTER ERMEAN Enter Mean of Error 
‘ ERRDEV Deviation of Error 
> \ 
[i RTMEAN Mean Error Rate 
[i RTDEV Deviation Error Rate 
! f FORMAT Whenever a format is called for, 


parentheses must be used. 


